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ABSTRACT 


Two  examples  of  quantum  mechanical  dissipative 
systems  are  studied.  The  first  one  is  the  multichannel 
Schrodinger  equation  with  a  nonlocal  separable  potential. 

This  problem  can  be  solved  exactly  and  the  partial-wave  scatter¬ 
ing  amplitude  can  be  found.  A  comparison  is  made  with  the 
well-known  N/D  method  for  one-channel  and  for  multichannel 
problems.  It  is  found  that  the  generalization  of  treating 
N  and  D  as  matrices  in  the  n-channel  problem  is  not 
possible  for  this  solvable  example.  Bound-state  and  resonance 
poles  of  the  partial-wave  amplitude  are  also  discussed. 

The  second  example  is  that  of  a  bound  electron 
(assumed  to  have  two  levels  only)  interacting  with  the 
electromagnetic  field.  This  model  is  similar  to  the  Lee 
model.  In  this  model,  the  transition  of  the  bound  electron 
from  the  higher  level  to  the  lower  level  is  accompanied  by 
the  emission  of  one  photon.  Results  for  level  shift,  line 
breadth  and  the  probability  amplitude  are  obtained,  which 
correspond  to  the  characteristics  of  a  classical  damped 
system,  i.e.  natural  line  breadth  and  the  exponential  decay 
form.  Modification  to  this  model  is  also  considered. 

This  modification  essentially  changes  the  eigenstates  of 
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the  photons,  i.e.  the  bare  photon  becomes  a  coherent  mixture 
of  physical  photons.  Hence,  the  original  interaction  term 


is  altered  accordingly. 
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Chapter  1. 


INTRODUCTION 


In  classical  mechanics,  when  the  forces  acting  on 
a  particle  (or  particles)  are  known  we  can  write  down  the 
equation  of  motion  by  virtue  of  Newton's  Second  Law  of  Motion. 
With  the  given  initial  conditions  on  position  and  velocity  of 
the  particle,  the  equation  of  motion  can  be  integrated  and  thus 
the  motion  of  the  particle  is  determined  completely.  We  note 
here  that  the  force  on  the  particle  can  be  a  function  of 
position  and  velocity  of  the  particle  as  well  as  an  explicit 
function  of  time,  i.e.  the  force  may  either  be  conservative 
or  dissipative. 

One  can  introduce  canonical  variables  to  form  a 

Hamiltonian  function  which  is  the  total  energy  of  the  particle 

(expressed  in  terms  of  canonical  variables).  With  the 

Hamiltonian  function,  the  classical  equation  of  motion  can 

1 

be  written  concisely  in  Poisson  brackets  .  A  natural 
way  to  quantize  the  classical  equation  of  motion  is  to  treat 
the  canonical  variables  as  dynamical  operators  (usually 
noncommutat ive )  in  quantum  mechanics  and  to  interpret  the 
Poisson  bracket  as  l/ifi  times  the  commutator.  In  this 
way,  the  classical  equation  of  motion  can  be  taken  into 
quantum  theory  by  Heisenberg's  equation  of  motion.  Another 
way  is  by  the  use  of  Schrodinger  equation  whose  solution 
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corresponds  to  the  Hamilton's  principal  function  in  the 
geometrical  optics  limit.  These  two  methods  are  mathematically 
equivalent2 


When  one  looks  at  the  Heisenberg  picture  of 

quantization,  the  imposition  of  the  Poisson  bracket  on  the 

commutator  would  seem  an  added  postulate  in  the  theory. 

This  raises  the  question:  "Do  the  equations  of  motion 

determine  the  quantum  mechanical  commutation  relations?" 

(by  Wigner)  .  The  answer  is  found  to  depend  on  the  form  of 

the  Hamiltonian  and  is  in  the  negative  for  a  free  particle 

and  for  the  simple  harmonic  oscillator.  This  nonuniqueness 

in  the  commutation  relation  for  the  oscillator  has  been 

4 

resolved  by  several  authors  .  Granting  the  Heisenberg's 
equations  of  motion  for  operators  determine  the  commutation 
relations  uniquely,  there  is  still  the  question  of  the  form 
of  Hamiltonian. 

Usually,  the  total  energy  of  the  system,  when  expressed 
in  canonical  variables,  is  the  Hamiltonian  function  which  can 
then  be  taken  over  as  the  Hamiltonian  operator  in  Heisenberg 

5 

picture.  However,  as  pointed  out  by  Havas  ,  for  a  given 
classical  equation  of  motion  many  different  Hamiltonians 
exist,  one  of  which  may  be  the  total  energy  of  the  system. 
Quantization  procedures  using  different  Hamiltonian  other 
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than  the  total  energy  will  be  quite  different  from  the 
Heisenberg  method  or  impossible,  even  though  all  these 
Hamiltonians  do  generate  the  same  classical  equation  of 
motion.  At  this  point,  one  would  rather  take  the  notion 
that  equations  of  mot ion, classical  or  quantal,  are 
fundamental  and  unique,  then  one  would  have  to  look  for 
different  equations  for  different  systems.  This  can  be  done 
in  classical  mechanics,  but  one  is  lost,  when  dealing  with 
quantum  mechanics,  in  searching  for  the  equation  of  motion 
without  falling  back  on  the  Hamiltonian. 

In  this  thesis,  we  are  interested  in  dissipative 
system.  For  instance,  in  classical  dynamics,  the  simple 
equation  of  motion  for  a  damped  oscillator  can  be  solved 
uniquely  with  two  given  initial  conditions  on  position  and 
velocity  of  the  oscillator.  When  we  try  to  study  this 
problem  in  quantum  mechanics,  we  have  to  find  the  Hamiltonian 
that  generates  Heisenberg's  equations  of  motion.  However, 
since  this  Hamiltonian  is  not  the  total  energy,  the  usual 
Heisenberg's  quantization  procedure  does  not  apply  and 
the  correct  quantization  rule  for  this  particular  Hamiltonian 
is  not  known.  Here,  Hamilton's  canonical  equations  are 
inconsistent  with  the  usual  commutation  relations.  Another 
way  to  consider  this  problem  is  just  to  "write  down"  the 
correct  quantal  equation  of  motion  for  the  damped  oscillator 
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avoiding  the  use  of  Hamiltonian  entirely.  Of  course,,  we 
search  in  vain  for  this  quantal  equation.  There  seems  to  be 
a  third  way  of  trying  to  solve  this  problem  by  taking  the 
total  energy  as  the  Hamiltonian  (which  is  not  explicitly 

c 

time-dependent)  plus  the  Rayleigh's  dissipation  function 

This  last  hope  diminishes  quickly,  as  pointed  out  by 
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Brittin  that,  for  a  dissipative  system  with  a  Hamiltonian 
not  explicitly  depending  on  time,  quantization  can  be 
carried  out  in  Heisenberg  picture  only  if  the  dissipative 
force  is  a  function  of  position  and  no  dissipative  force 
is  possible  in  Schrodinger  picture.  Since  we  believe  in 
the  equivalence  of  these  two  pictures  (two  different  methods 
of  quantization),  the  above  inconsistency  rules  out  the  hope 
to  quantize  a  classical  system  with  velocity-dependent 
dissipative  force  like  that  of  the  damped  oscillator 
(classical  electrodynamics  is  the  exception).  We  accept 
the  idea  that  a  dissipative  system  can  only  be  treated  as 
a  subsystem  of  a  conservative  system  in  quantum  mechanics, 
in  contrast  to  classical  mechanics  in  which  both  conservative 
and  dissipative  systems  can  be  studied  separately  from 
their  respective  equations  of  motion. 

We  consider  some  examples  of  dissipative  systems 
in  quantum  mechanics.  The  multichannel  problem  (Chapter  3.) 
has  been  applied  to  inelastic  nuclear  reactions  for  some 
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years  .  The  effect  of  all  other  channels  on  a  particular 

channel  is  to  create  an  absorptive  potential  which  can  be 
considered  as  the  quantum  mechanical  analogue  of  a  damped 
classical  system.  We  choose  a  solvable  example  with  nonlocal 
separable  potential  to  study  mainly  the  behaviour  of  the 
S  matrix  and  to  compare  the  exact  form  of  the  S  matrix 
with  the  well-known  N/D  method  used  in  dispersion  relations. 


Another  example  (Chapter  4.)  is  that  of  a  bound  charged 

particle  interacting  with  an  electromagnetic  field  (non- 

relat ivistic ) .  A  bound  electron  (having  two  levels  only) 

decays  into  the  lower  level  with  the  emission  of  one  photon 

and  the  reverse  process.  This  problem,  which  in  some  ways 

10 

is  similar  to  the  Lee  model  ,  can  be  solved  exactly.  We 
obtain  results  for  level  shift  and  line  breadth  which  can 
be  compared  with  those  obtained  by  time-dependent  perturbation 
method.  More  important  is  that,  in  classical  electrodynamics, 
a  moving  point  charge  emits  radiation  which  in  turn  reacts 
on  the  charge's  own  motion.  The  equation  of  motion  for  the 
charge  (Dirac's  equation)  is  in  certain  ways  similar  to 
that  of  a  damped  harmonic  oscillator.  The  emitted  line 
by  the  charge  is  not  infinitely  sharp  but  has  a  certain 
natural  breadth  due  to  the  damping  force  of  the  emitted 
radiation  on  the  charge  itself  (Heitler's  §  4.)^  .  Now 

our  quantum  mechanical  dissipative  system  (which  is  a  sub- 
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system  of  a  conservative  system)  exhibits  the  same  characteris¬ 
tics  a  classical  damped  system.  This  convinces  us  that 
although  a  classical  damped  system  cannot  be  quantized 
directly,,  it  can  always  be  considered  as  part  of  a  conservative 
system  which  can  be  quantized,  and  then  the  quantum  mechanical 
damped  behaviour  (quite  similar  to  that  of  classical)  can  be 
found  from  the  complete  solution  of  the  conservative  system. 

Modification  to  the  above  model  is  presented  in 
Chapter  which  changes  the  interaction  term  as  well  as 
the  vacuum  state  for  the  photon. 


sd 
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Chapter  2.  QUANTIZATION  OF  DAMFED  CLASSICAL  SYSTEM 

A  damped  classical  system  can  be  specified  uniquely 
by  its  equation  of  motion  (plus  initial  conditions).  Generally, 
the  classical  damped  system  cannot  be  quantized  directly.  In 
Section  1,  we  discuss  some  of  the  difficulties  in  quantizing 
a  damped  classical  system  and  conclude  that  quantum  mechani¬ 
cally  a  damped  system  can  only  be  treated  as  a  subsystem  of 
a  conservative  system.  We  present  a  simple  quantal  system 
(Section  2.),  similar  to  the  Wigner-Weisskopf  model,  which 
exhibits  the  character  of  damping. 

1 .  Quantization  of  Damped  Classical  System 

5 

It  has  been  shown  by  Havas  that  with  a  suitable 
"integrating  factor"  a  very  broad  class  of  classical  equations 
of  motion  can  be  derived  from  a  variational  principle  and 
hence  the  Lagrangian  and  the  Hamiltonian  can  also  be  obtained 
in  the  usual  way.  However,  due  to  the  multiplicity  of  possible 
integrating  factors,  for  a  given  equation  of  motion  many 
different  Lagrangians  exist.  Quantization  for  different  Lagran¬ 
gian  will  be  quite  different  or  even  impossible.  In  spite  of 

12 

this  general  finding,  there  are  several  papers  attempting 
to  quantize  the  damped  harmonic  oscillator.  The  common 
feature  in  those  papers  is  that  starting  from  the  classical 
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equation  of  motion  for  the  damped  oscillator, 

x  +  ax  +  oo2x  =  0,  (2.1) 

which  is  derivable  from  a  time-dependent  Hamiltonian  (not  the 
energy  of  the  damped  oscillator), 

H  =  (l/2m)e  atp2  +  hmu2ea^x2 ,  (2.2) 

where  p  me<xt  x.  Quantization  is  then  by  the  usual  commu¬ 
tation  relation 


[x5  p]=  ih. 


(2.3) 


or 


[x,  x]  =  fhe  at/m. 


(2.4) 


Eq.  (2.4)  implies  that 

AxAx  >(li/m)e-a1:  (2.5) 

violates  the  uncertainty  principle  for  an  oscillator. 


AxAx  >  Jfi/m , 


(2.6) 


. 
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which  is  valid  even  for  the  damped  oscillator 

The  violation  of  the  uncertainty  principle  may 

seem  to  arise  from  the  explicit  dependence  of  H  and  p  on 

time.  Some  years  ago  Seeger  considered  the  same  problem. 

He  chose  a  time-independent  Hamiltonian  plus  the  Rayleigh’s 

dissipation  function  (depending  on  velocity)  and  then 

quantized  the  damped  oscillator  by  the  ordinary  commutation 

relation  working  entirely  in  the  matrix  form  of  quantum 

mechanics,  Heisenberg  representation.  However,  it  was  shown 

7 

by  Brittin  that  for  a  dissipative  system  not  explicitly 

depending  on  time,  quantization  can  be  carried  out  in  the 

Heisenberg  representation  only  if  the  dissipative  force  is 

a  function  of  position  (not  velocity)  and  no  dissipative 

force  which  is  a  function  of  position  or  momentum  is  possible 

in  Schrodinger  representation.  This  kind  of  inconsistency 

rules  out  the  possibility  of  quantizing  a  dissipative  system 

by  the  usual  procedure.  Even  within  the  frame  work  of  quantum 

13 

mechanics,  Razavy  recently  showed  that  in  nuclear  physics 
an  often  used  class  of  velocity-dependent  Hamiltonians, 
which  attempts  to  explain  the  strong  short-range  repulsion 
of  the  two-nucleon  interaction  at  high  energies,  leads  to 
unacceptable  results.  Though  the  velocity-dependent 
Hamiltonian  satisfies  the  general  requirements  of  Hermiticity 
and  invariance  properties,  the  Heisenberg  equation  of  motion 
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for  the  energy  operator  (closely  related  to  the  Hamiltonian) 
cannot  be  satisfied. 

In  view  of  the  difficulties  of  quantization,, 

apparently  the  Hamiltonian  chosen  for  quantization  can  only 

5 

be  the  energy  of  the  system  which  is  uniquely  defined  in 
classical  mechanics.  As  for  the  quantization  of  a  dissipative 
system,  one  should  treat  the  dissipative  system  as  a  subsystem 
of  a  non-dissipative  system.  The  complete  system  will  then 
be  solved  by  some  approximation  if  needed.  The  behaviour 
of  the  dissipative  system  can  then  be  interpreted  from  the 
solution  of  the  complete  non-dissipative  system. 

2 .  The  Wigner-Weisskopf  Model 

A  simple  solvable  model  for  inelastic  processes 

14 

was  considered  some  years  ago  by  Wigner  and  Weisskopf 
It  consists  of  a  motionless  particle  whose  wave  function  is 
X(t),  and  another  moving  particle  of  mass  m  with  the  wave 
function  Y/(r,  t).  These  two  particles  act  as  one  another's 
sources  with  a  real  form  factor  p(r).  The  Schrodinger 
equations  are 


(i8/3t  +  V2/2m)^(r ,t )  =  Gp(r)x(t) 


^Ino  nso  noltfssitasjjp  nssorio  rreino^IiaisH  er  y,  '  91 
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(id/dt  -  li)x(t) 


r_ 

G  p(r)^(r,t)dr. 


(2.7) 


Assuming  p(r)  depends  on  the  magnitude  of  the  radius 
vector,  then  only  S-waves  are  coupled.  Usually  p( r)  is 
of  a  very  short  range  rn  with  the  point  interaction  (local) 

limit  as  rQ  0.  The  solutions  for  (2.7)  have  terms  in  r0_1  . 

Vb 

Martin  studied  the  time  development  of  a  prepared  state 
|  ip  =  0,"X_  =  1^  for  this  model  at  t  =  0.  He  found  that  such 
a  state  is  not  physically  acceptable,  since  the  mean  energy 
of  the  system  diverges  at  t  =  0,  which  implies  that  the 
decay  curve  has  a  cusp  at  t  —  0.  Recently  Razavy16  remedied 
the  above  defect  by  a  noncausal  coupling  in  the  Schrodinger 
equations.  In  the  following  we  shall  consider  a  similar 
system  with  nonlocal  interaction,  which  will  give  results 
similar  to  Razavy’ s. 

The  Schrodinger  equations  for  two  channels  coupled 
by  a  nonlocal  separable  potential  are 

(V2+k2)^(r)  =  GV(r)  |v(r’  )<f>g(r  '  )dr  ’ 

(  V2  -k2  )  (p  (r)  =  G  V  ( r )  (v  (  r  ’  )  b  ( r  ’  )  dr  ’  , 

i  J  k 


(2.8) 


. 


i  ' 

. 
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where  |C2  =  £  -  k2"  >  0,  £  is  the  internal  energy  of  the  two 
particles  described  by  the  wave  function  <j>  .  Eq.  (2.8) 
will  be  considered  in  more  detail  in  the  next  chapter.  Our 
main  purpose  here  is  to  show  that  the  phase  shift  in  ^  is 
qualitatively  the  same  as  Razavy's  and  a  state  of  finite 
mean  energy  can  be  prepared  for  the  system  at  t  =•  0. 

Since  V(r)  depends  on  |r|  ,  Eq .  (2.8)  will  only 
couple  S-waves.  We  have  assumed  that  K2  >  0,  so  there 
will  be  no  incoming  wave  in  <f>  as  well  as  no  spherical 
outgoing  wave,  i.e.  <j>  — »  e  /r  as  r-^co.  With  little 
modification,  namely  by  putting 

gn=g22-4  gi2  =  ^irG,  k?=k2,  and  kl  =  -K2  ,  (2.9) 

the  results  for  the  wave  functions  and  the  partial-wave 
amplitude  in  chapter  Eq.  (3«ll)j  (3.12),  and  (3 - 13 ) ^  can 

be  taken  over  to  give 


p2-k2-ie  j 


V2(p,)dp»/, 


p  ’  2  + 


K‘ 


r  r 

■ 
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<M?) 


G  V  ( k )  V  ( p  )  ( 

P2  +  K2  / 


(2.10) 


where 


A  =  1-G: 


J 


vt  p ) dp 


p2-k2-ie 


vt  p  1  )  dp  '  . 


(2.11) 


7 


P  ’  2  +K  2 


The  tilde  indicates  Fourier  transforms  of  the  respective 

quantities.  We  may  choose  V(p )  l/(kz  +■  b2"),  Eq  .  (3.18), 

✓ — " 

then  from  ^(p),  we  can  find  ^^.(r)  whose  asymptotic 

k  k 

condition  at  r-*oo  yields  the  scattering  amplitude  (  S-wave), 


a 


2tt 2G 2 V2  ( k ) 

A 


J 


V2 (p )dp 
P2+k2 


k 


(2.12) 


' 
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From  this  relation  it  follows  that 


kcotS  =  Re( a  1 ) 


r 


[1-G 


2 

> 


V2 (p ' )dp * 

p ’ 2  +  k2 


X 


[2tt2G2V2  (k) 


V2 (p)dpj -1 

P2  +  k2 


(2.13) 


where  P  stands  for  the  principal  value  of  the  integral. 

Eq.  (2.13)  is  the  ratio  of  two  polynomials  in  k2  and 

2_ 

Razavy's  result  is  a  polynomial  in  k  .  Hence,  our  system 
gives  qualitatively  the  same  scattering  as  the  one  with 
noncausal  coupling.  This  is  to  be  expected.  As  pointed  out 
by  Gasiorowicz  and  Ruderman^  that,  under  certain  conditions, 
nonlocal  interaction  is  equivalent  to  noncausal  interaction, 
in  the  sense  that  they  produce  the  same  scattering  amplitude. 

For  simplicity,  we  shall  assume  there  is  only  one 

bound  state  4>  (normalized)  in  <p  when  G  -  0.  Let  us 

1  b 

prepare  a  state  |''Jf(t  =  0)^>  =  1 1^=  0,  <p  y  for  our  system 

at  t  =•  0  and  study  its  subsequent  development.  We  can 


;r 

, 
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expand 


with 


li(o)>  in  terms  of  the  eigenstates  (2.10), 


|y(o)> 


r 


dkX(k)  |  r )  3 4>_^ ( r* ) >  , 
k 


■+ 

k 


(2. 14) 


Mk)  given  by 


X(k)  =  <\p  (?)  ,<j>  (?)  1 0,4>  > 
k  k  u 


r 


=  (  2tt) 


-3/2 


<j>^(p) 

k 


IP 


dp(J)b(r)dr 


r 


=  (  2  7T  ) 


•3/2 


) 


.  ->■  —  K  x1 

p  )  e^p  ‘  rdp(  /Kb/2TTe  13  /r)dr 
k 


=  /k^/tt \<j>  (p)/(p: 


+ 


K^)dp 


(2.15) 
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At  a  later  time,  the  state  will  be  given  by 


nt)> 


f 


->  . 


dkX (k) e 


-ik2 1 


K?)  , 

l 


6  (r)>  . 

I 


(2.16) 


The  mean  energy  of  the  system  at  t  =  0  is 


< 


:¥(  0)  | H 


no) 


where  H  is  the  Hamiltonian  (a  matrix)  for  the  two-channel 
Schrodinger  equations  (2.8).  With  /V(k)  given  by  (2.15), 
the  mean  energy  above  is  finite.  Hence,  the  decay  curve  for 
this  model  does  not  have  a  cusp  at  t  =■  0.  The  probability 
amplitude  that  the  system  remains  in  1^(0))  after  t  is 
determined  by  ■('y(0)|  T|f(t)\  which  has  damped  exponential 
behaviour  for  small  t  and  ultimately  goes  as  some  inverse 
power  of  t  for  t  -^°o  .  More  detail  of  a  similar  prepared 
state  will  be  given  in  Chapter  4. 


, 
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Chapter  3-  MULTICHANNEL  PROBLEM  WITH  NONLOCAL 

SEPARABLE  POTENTIAL 

A  well-known  method  of  treating  inelastic  scattering 
and  resonance  reaction  is  to  extend  the  one-channel  Schrodinger 
equation  to  the  multichannel  equation.  After  a  brief  discussion 
on  the  multichanned  Schrodinger  equation  (Section  1.),  we 
solve  an  example  with  a  nonlocal  separable  potential 
(Section  2.).  Although,  for  one  channel,  the  partial-wave 
amplitude  can  be  put  into  the  often  used  N/D  form,  a  straight¬ 
forward  generalization  of  N/D  to  multichannel  is  not 
possible  for  this  solvable  example  (Section  3«)*  In  Section 
4,  the  analytic  properties  of  the  2-channel  partial-wave 
amplitude  is  studied  with  the  attention  to  bound  states  and 
resonances . 

1 .  The  Coupled  Schrodinger  Equation 

A  general  situation  in  nuclear  scattering  experiment 

consists  of  a  projectile  nucleon  impinging  on  a  target  nucleus 

which  can  be  excited  into  many  higher  levels.  Such  a  reaction 

8  9 

can  be  described  by  the  coupled-channel  Schrodinger  equation  5  : 

-^M_1V2?  +  V?  = 


(3-D 
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where  we  consider  n  coupled  channels  interacting  with  two- 
body  potentials,  Eq.  (3.l)  is  written  in  the  centre-of -mass 
coordinates  of  the  two-body  channels.  The  name  "channel" 
has  the  following  meaning:  the  incoming  nucleon  plus  the 
target  nucleus  in  its  ground  state  is  designated  the  incident 
channel  which  is  capable  of  going  into  n  exit  channels 
having  an  out-going  nucleon  plus  the  target  nucleus  in  one 
of  its  excited  states  or  its  ground  state.  We  confine  our¬ 
selves  to  two-body  interactions  by  V  so  that  particles  in 
each  channel  are  grouped  into  two  bodies,  e.g.  projectile 
plus  target,  out-going  particle  plus  excited  nucleus.  The 
diagonal  matrix  M  has  elements  ><*,,  the  reduced  mass  of 
the  two  bodies  in  a  certain  channel  <*,.  The  matrix  £,  has 
element  s  -  £<*,  )  y  where  E  is  the  total  energy 

of  the  whole  system  or  the  kinetic  energy  of  the  incident  nucleon 
(taking  ground  state  energy  of  the  target  as  zero),  and  £oc 
is  the  internal  energy  of  the  two  bodies  in  channel  cL  or 
the  excitation  energy  of  the  target  nucleus.  We  shall  call 
Ex.  the  "threshold"  of  channel  PC  because  it  is  the 
minimum  kinetic  energy  that  the  incident  nucleon  must  possess 
in  order  to  make  the  interaction  energetically  possible  in 

channel  oC  .  The  potential  matrix  V  is  symmetric  due  to 

18  — 

time-reversal  invariance  .  Y  is  the  wave-function 
matrix  describing  the  interaction  between  n  channels.  Each 
column  of  ^  and  the  whole  matrix  are  solutions  of 


,  . 


-I  ' 
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(3.1). 


Eq.  (3.l)  can  be  reduced  to  a  simpler  form  by 

defining 


V 


2M2  V  M' 


(3.2) 


±  IA>  +  I/O 

where  M  is  the  diagonal  matrix  with  elements 
Substitution  of  (3.2)  into  (3.1)  yields 

-V2Y  +  VT  =  K2^  (3-3) 

where  K  is  the  diagonal  matrix  of  the  channel  wave  numbers, 
—  / 2  (E  -  g^).  We  may,  therefore,  consider  all 

channels  having  the  same  reduced  mass  1/2,  and  take  (3-3) 
as  our  coupled-channel  equation.  At  a  given  total  energy  E 
it  may  be  possible  to  excite  the  target  nucleus  into  some 
levels  and  still  have  some  (positive)  kinetic  energy  left 
over  for  the  out-going  particle;  other  higher  levels  may  not 
be  so  accessible.  The  channels  corresponding  to  the  former 
are  said  to  be  open,  i.e.  E  >  g^,  and  those  corresponding 
to  the  latter  are  said  to  be  closed,  E  <  g^ .  The  energy 


■■  '  ■ 


, 

.  ' 
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at  which  a  channel  opens  up  is  its  threshold. 

2 o  A  Solvable  Example  of  Nonlocal  Separable  Potential 

19 

It  is  well  known  that  effects  of  many-body  forces 
exchange  forces  (Majorana  ,  particle-exchange  ),  velocity 
dependent  forces,  and  the  usual  relativistic  correction  to 
the  energy,  all  give  rise  to  nonlocal  potentials  in  the 
Schrodinger  equation.  In  the  coupled  equation  (3.3)  with 
local  (real)  interaction,  we  can  single  out  a  particular 
channel,  e.g.  the  incident  channel  1,  by  eliminating  all 
other  channels.  The  resulting  Schrodinger  equation  for  chan¬ 
nel  1  contains  a  generalized  "optical"  potential^  whose 
imaginary  part  is  negative  definite.  This  means  that  the 
effect  of  other  channels  on  1  is  to  produce  an  absorptive 
potential  since  the  incident  nucleon  may  be  "absorbed" 

(not  literally  but  with  kinetic  energy  changed)  in  different 
ways  through  n  exit  channels.  This  generalized  "optical" 
potential  is  nonlocal  though  not  separable  *.  This  is  the 
quantum  mechanical  analogue  of  the  damped  classical  system. 

In  our  discussion  we  shall  take  a  real  nonlocal 
separable  potential  matrix,  which  will  be  specified  later. 


*  A  general  type  of  nonlocal  potentials,  with  separable 
and  local  as  two  extreme  cases,  can  be  expressed  as  a  sum  of 
separable  potentials.212' 


■' 


:.1 


■ 
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and  shall  ignore  spin  variables  in  the  coupled  Schrodinger 
equation  (3. 3).  Let  us  make  the  usual  partial-wave  decomposi¬ 
tion  for  the  wave  function  and  the  nonlocal  potential. 


=  £uie)(r)/r  Vcose) 

X/ 


V  (?,?')  =  l  (  2£+l)V^  (r  ,r  1  )/(  ^Trrr  '  )  P0(f*f')  ,  (  3 .  *0 

ay  0  a  y  x, 


where  we  have  assumed  the  nonlocal  potential  to  be  "central", 

->  — * 

i.e.  depending  only  on  the  magnitude  of  the  distance  r-r' 

(r  and  r1  are  the  unit  vectors  in  the  directions  of  r  and  r' 
respectively).  The  nonlocal  potential  is  not  only  symmetric 

in  oL  and  y  ,  as  mentioned  in  the  last  section,  but  also 

_  23 

symmetric  in  r  and  r'  due  to  the  requirement  of  Hermiticity. 

Then  Eq.  (3.3)  becomes 


[d2/dr2  +  k2  -  £  (  £+1  )/r2  ]  ^  (r) 

o t  a\5 


n 


- 1 
Y=1 


r 


O 


v(£)  (r  r  ’  (r!  )dr  ’ 

ay  s  y£ 


0 


(3.5) 


. 


.  r  ' 

. 


. 
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which  shows  that  only  partial  waves  of  the  same  order  are 
coupled.  For  our  real  nonlocal  separable  potential,  we  choose 
the  simple  form, 


„(£)/  *  \  (£)  (£),  v  U),  •  \ 

(r5r')  =  v  (r ) v  '(r'), 


(3.6) 


where  -  g  are  the  real  coupling  constants.  The 

OL  7  f  cL 

S  matrix  for  the  1-th  partial  wave  is  determined  by  the 


boundary  condition  at  r  oo 


24 


e-^-e^^e-i^S^b.  (3.7) 

a3  a  3  a$  ap 


In  the  following  discussion  we  shall  suppress  the  JL  indices 
(1) 


in  u^p  ,  g^|  ,  v(£)  ,  and  S 


(I) 

* 


Eq.  ( 3  -  5 )  with  the  potential  given  by  (3.6),  can 
be  written  in  momentum  space  by  the  following  transform: 


( 


o° 


u  (r)/r  =  /2/tt 
a  p 


Ua3(p)j*  £(pr)p2dp 


J 


OO 


r 


v(r)/r  =  /2/ 


TT 


J 


v(p) j  £(pr)p2dp 


. 

I  1 


f! 
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CO 


uaB(p)  =  1/577 


J 


ua^(r)/r  j£(pr)r2dr 


O© 


v ( p  )  =  /2/tt 


v(r)/r  j  ^ (pr )r2dr  . 


(3.8) 


From  (3.6)  and  (3.8),  we  obtain  for  ( 3 - 5 ) 


oO 


(ka"P  )uae(p)-2SaYV(p) 


7 


v(p ’ )u  g(p ’ )p  f  2 dp  '  =  0 


(3-9) 


with  the  general  solution: 


6(k  -p) 
a 


uaB(p) 


=  6 


a£ 


v(p) 


r 


OO 


k2 

a 


p2-k2-ie 
^  a 


•Is 


ay 


v(p ’ )u^^(p ’ )p ’ 2dp '  .  (3.10) 


The  integration  constants  J  v(p)  Uyg(p)p2-dp  can  be  evaluated 

o  ' 

from  a  set  of  simultaneous  equations  obtained  by  substitutions 
of  (3.10)  into  the  constants  themselves*  We  get  finally 


u  „(p)  =  6  aS.(1^-P) - 2^- 

aB  aS  k2  p  2  — k  2  — i £ 

a  ^  a 


Tg  (A  1)  Qv(k  )  (3.11) 

2  lr2  A  ay  y3  ah 
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where  the  matrix  A  is  given  by 


A 


a3 


oo 

r 


+ 


v2 (p)p2dp 

p2-k2-ie 

a 


(3.12) 


The  S  matrix  can  be  obtained  by  writing  u^^p) 
in  configuration  space,  Eq.  (3.8),,  and  from  the  asymptotic 
condition  ( 3 - 7 ) 


a3 


=  6 


a3 


iir/k  k  v(k  )v(k  )7g  (A  1)  D 
a  3  a  3  L  ay  y3 


=  6 


a3 


2i/k  kD 
a  3 


a3 


(3.13) 


where  a^  is  the  i -  th  partial-wave  amplitude.  As  implied 

by  time-reversal  invariance,  the  S  matrix  is  symmetric, 

which  can  also  be  shown  explicitly  from  (3*13)  and  (3.12). 

We  note  that,  from  (3.12),  A^  can  be  considered  as  a 

2. 

function  of  the  total  energy  E  =  k^  -f  ,  which  may  also  be 

analytically  continued  for  a  complex  E  in  a  certain  region 
if  we  choose  v'2  (p)  analytic  in  that  region  (and  of  course 
vanishing  fast  enough  for  large  p  to  make  the  integral  in 
(3.12)  convergent).  As  a  function  of  the  complex  variable  E, 


. 


t  .  ..  . 

A  '  '  , 
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A^a  will  have  a  cut  on  the  positive  real  axis  starting  from 
the  threshold  £<*,  to  infinity  (the  ground  state  energy 
can  be  taken  as  zero).  Since  we  have  taken  vr2~  (p)  to  vanish 
for  large  p,  A^— >  <5^  as  E  -^<*>  ,  i.e.  A ^  is  analytic  and 
normalized  for  E  —  oo  .  in  the  following  sections,  we  shall 
study  the  analytic  properties  of  the  S  matrix  in  more  detail. 

3 •  The  N/P  Form  for  the  Partial-Wave  Amplitude 

It  has  been  shown  by  many  authors2^  that,  for  a 

one-channel  problem,  the  partial-wave  amplitude  a(s)  is  a 

function  of  s,  the  total  centre-of-mass  energy  squared. 

It  has  a  right-hand  cut  (kinematical  cut)  in  the  complex  s 

plane  coming  from  the  unitarity  of  the  S  matrix  and  a  left- 

hand  cut  (dynamical  cut)  due  to  the  interaction  between  the 

two  particles  in  the  channel.  For  nonrelativist ic  potential 

scattering  D  that  we  are  interested  in,  s  =  k  the  centre- 

of-mass  momentum  squared  if  we  take  the  internal  energy  of 

the  two  particles  as  zero,  i.e.  the  threshold  for  this 

2_ 

interaction  is  at  k  =  0.  The  general  analytic  behaviour 

mentioned  above  enables  one  to  write 


j 


a  =  N/D 


(3.14) 


. 

i 
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where  D  has  the  same  right-hand  cut  as  a  but  no  left-hand 
cut,  while  N  has  the  same  left-hand  cut  as  a  with  no 

right-hand  cut.  D  is  normalized  at  k2  =  00  and  its  zeros 

with  Im(k)  >  0  (Re  k  =  0)  correspond  to  the  bound  states. 

N  and  D  can  then  be  expressed  as  some  kind  of  dispersion 

relations.  These  two  dispersion  relations,  two  coupled 
integral  equations,  may  be  solved  by  successive  approximations 
This  is  a  roundabout  way  to  solve  a  scattering  problem,  though 
we  do  not  need  it  to  solve  potential  scattering  since  it  is 
easier  to  solve  the  Schrodinger  equation.  However,  at  high 
energies,  where  the  governing  equations  are  not  solvable  or 
simply  are  nonexistent,  the  N/D  method  furnishes  an  appro¬ 
ximation  that  can  be  handled  without  knowing  the  exact  dynam¬ 
ical  structure  of  the  system.  Usually  when  the  projectile 
gets  more  and  more  energetic  it  is  capable  of  producing  diff¬ 
erent  kinds  of  reactions,  elastic,  inelastic,  rearrangements. 
In  these  cases,  a  single-channel  formulation  will  not  be 
sufficient.  One  must  take  account  of  the  channels  that  are 
most  relevant  to  the  problem  concerned. 

There  have  been  many  attempts  to  generalize  the 

27  28  2 Q 

single-channel  N/D  method  to  multichannel  processes  5  5 

An  obvious  way  to  generalize  (3*1*0  is  to  treat  N  and  D 
2  7 

as  matrices  because  the  partial-wave  amplitude  becomes  an 

n  x  n  matrix  a  n  for  n  coupled  channels: 

a3 

aaB  =  ^  Ncr/D 


(3.15) 


' 
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with  the  usual  respective  analytic  properties  assigned  to 
the  N  and  D  matrices,  namely  that: 

(a)  dynamical  cut  only  in  N; 

(b)  physical  cut  (unitarity  cut)  only  in  D;  and 

(c)  bound-state  poles  are  produced  by  zeros  of  D  (the 
determinant  of  matrix  D  ) . 

Another  way  to  account  for  the  effect  of  other  channels  on  a 

particular  channel  is  to  introduce  an  "absorptive  parameter" 

29 

into  the  single-channel  scattering  amplitude  and  thus 

keep  the  simple  form  of  single-channel  N/D.  However,  it 

has  been  shown  that  results  obtained  from  "absorptive 

parameter"  N/D  are  not  equivalent  to  those  from  multichannel 
/  30 

N/D  „  Though  D  still  has  only  the  right-hand  cut,  N  will 

have  both  the  left  and  the  right-hand  cut.  Moreover,  condi- 

71 

tion  (c)  is  not  fulfilled 

For  our  solvable  example,  it  is  easy  to  check  if  the 
N/D  form  (single  or  multichannel)  is  acceptable  with  the 
above  three  conditions  imposed  .  For  one  channel,  from 

(3-13)  and  (3.12),  we  obtain  the  partial-wave  amplitude, 

a(k2 )  =  N(k2)/D(k2) 


oo 


r- 


=  — ^TTgV2  (k)/(l+g 


J 


v  2  ( P  )  P  2  dp 

p2-k2-ie 


) 


(3.16) 


' 

,  •  -  '  ;  1  ,  :  .  •  I  -  . 


. 


v  .  ! 


_ 
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which  fulfills  conditions  (a),, 


(b),  and  (c). 


provided 


32 

that  J 


v(r )/r 


0(r 


for  r  -*  0 , 


and 


v(r )/r~const .  exp(-br)  for  r  ->  °°  .  (3.17) 

We  may  then  assume  that  our  potential  (3.8)  is  of  the  form 

v(k)  -  l/(k2+b2),  (3.18) 

where  l/b  can  be  considered  as  the  range  of  the  potential. 
With  Eq.  (3.18),  the  left-hand  cut  of  N  becomes  an  isolated 
double  pole. 

For  multichannel,  we  have  already  noticed,  at  the 
end  of  last  section,  that  matrix  of  (3. 12)  (with  (3.17)) 

is  equivalent  to  the  D  matrix  in  (3.15).  If  we  make  the 
assumption  that  all  channels  have  the  same  threshold,  say  0, 
we  can  write  an  exact  N/D  form  for  our  partial-wave  amplitude, 
i .  e . 


aB 


=  l  N  (D  1)  Q 
L  ay  y8 


CO 

:  6  .  .  ■  '  -  -  .  -  .  1  •  •  ■ 


■  ■  ■  :  • 
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with 


(k) 


-W(k)gag 


o£> 

r 


D  Q(k)  =  A  Q(k)  =  6  Q- 
aB  aB  aB 


2 

TT 


Na6(P)p  dp 

p2-k2-ie 


(3-19) 


from  Eqs.  (3.12)  and  (3*13).  In  general, however,  different 
channels  will  have  different  thresholds,  as  the  target  nucleus 
would  have  excited  states  other  than  the  ground  state.  From 
(3*13)j  the  partial-wave  amplitude  (matrix)  can  be  written 
in  matrix  form 


a  =  -^TrvgA  1v  ,  (3-20) 

where  v  is  the  diagonal  matrix  with  elements  v^k^),  (contain¬ 
ing  the  left-hand  cut),  g  the  symmetric  matrix  of  the  coupling 
constants  gjQ  ,  and  matrix  A  is  given  by  (3.12).  Eq.  (3.20) 

00p 

is  no  longer  in  the  exact  N/D  form.  If  one  still  insists 
upon  writing  a^  =  2  Noty  (D~ '  )T/3  with  Dr/S  given  by  (3.12) 
a  different  N  matrix  will  be  found.  From  the  matrix 
equation  (3.20)  (a,  N,  and  D  are  now  matrices) 


. 
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a 


=  -^TTVgD 


N  = 


-^TTVgD 


(3.21) 


which  shows  that  N  will  have  the  right-hand  cut  (from  D) 
in  addition  to  the  usual  left-hand  cut  (from  "v ) .  This  may 
be  compared  with  the  situation  in  the  "absorptive  parameter" 
N/D22*^1  ,  except  that  the  zeros  of  det  (D)  in  (3.2l)  do 

produce  bound  states. 


In  conclusion,  a  straightforward  generalization  of 
the  N/D  form  for  the  partial-wave  amplitude  from  one-channel 
to  multichannel,  Eq .  (3.15)  plus  the  three  restrictions,  is 
not  applicable  to  our  simple  solvable  example,  and  we  doubt 
its  validity  in  the  much  more  complicated  problems  though  it 
has  already  been  used  in  the  treatment  of  strong  interactions. 

4 .  Analytic  Properties  of  the  Partial-Wave  Amplitude 

(  2  -  channel) 

We  shall  study  the  analytic  properties  of  the 
partial-wave  amplitude  in  the  simple  case  of  two  coupled 
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channels. *  The  general  behaviour  of  the  n-channel  amplitude 
is  similar  to  that  of  the  two-channel  amplitude.  The  matrix 
elements  a^  of  the  partial-wave  amplitude  can  be  found 
from  Eqs.  (3.13)  and  (3.12). 
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ai  i 


=  nn/A  = 


-7TV2  (ki  )  r  .  /  _  ~  2  \ 

Lg  1  1  +  V  g  1  1  g  2  2  — g  1  2  / 
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V2 (p)p2dp J 

p  2 -k2-ie 


a2  2  —  n22/A  - 


—  TT  V  2  (  k2  ) 

2  A 


[g22  +  (gl  1  g  2  2  — g  12)1 


v2 (p)p2dp 

p2 -k2 -ie 


] 


<v 

ai2  =  a2  1  =  n  1  2/A  =  -ttv( k  1  ) v ( k2  ) g  1  2/2A 


(3.22) 


where  A  is  the  determinant  of  the  matrix  A,  Eq .  (3. 12), 
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1  1 


J 
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PO 


oo 


r 


~gl  2 


v2 (p )p  2 dp 


p-kf-ie 


v2 (p ’ )p  ' 2 dp  * 

p ' 2 -k|-ie 


(3.23) 


Let  us  suppose  that  channel  2  has  a  higher  threshold  than 
channel  1,  i.e.  >£J(=0,  the  ground  state  energy). 

The  channel  momenta  ,  kx  and  k^  are  not  independent 
(their  signs  are  independent)  but  are  related  to  the  total 
energy  E  of  the  system  (the  kinetic  energy  of  the  incident 
particle)  by 


k?  =  E  ,  k|  =  E  -  e 2  j  (3.24) 

where  we  have  taken  the  channel  reduced  mass  as  1/2.  When 
aocp  are  considered  as  functions  of  the  complex  variable  E, 
they  have  (kinematical )  cuts  (from  the  integrals)  on  the 
positive  real  axis,  starting  from  the  thresholds  at  E  =  0 
and  E-  to  infinity,  and  (dynamical)  cuts  (from  the 
potential)  below  the  thresholds  along  the  real  axis  to  -oo 
We  also  note  that  a^-*  0*  and  ^  1*  for  |e|-*0°  .  In 
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the  following  discussion,  we  shall  follow  a  paper  by  Peierls 

The  characteristics  of  the  2-channel  processes 
depend  on  the  incident  energy  E  (real  and  positive).  If  E 


, 
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is  less  than f£ the  only  possible  real  process  is  elastic 
scattering  in  channel  1.  If  E  is  above  the  threshold 
for  channel  2,  inelastic  scattering  is  also  possible.  The 
zeros  of  A  ,  with  certain  restrictions  on  k^  and  k^  , 
corresponds  to  the  bound-state  and  resonant  poles  of  a^  . 
Their  positions  are  determined  from  the  (complex)  roots 
Of  (3.23). 


Let  us  look  at  the  bound-state  poles  first.  There 
are  two  kinds  of  bound  states,,  one  has  the  bound-state  energy 
below  all  the  thresholds  (the  ordinary  kind)  and  the  other 
has  the  bound-state  energy  between  the  two  thresholds.  The 
latter  kind  is  called  the  bound  state  embedded  in  the 
continuum^  because  it  looks  like  an  ordinary  bound  state 
in  channel  2  but  with  an  energy  in  the  continuum  of 
channel  1,  so  it  is  not  a  bound  state  in  channel  1.  Bound 
states  of  the  first  kind  are  given  by  the  roots  of  (3.23) 
with  both  kA  and  k2  positive  imaginary,  i.e. 
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(3-25) 
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where  k  t  =  i  (Ct  =  -t-i  /"e^  ,  k  ^  =  ilCz  =  +  i  -  <£2  ,  and 

E ^  is  the  negative  bound-state  energy.  It  is  seen  that 

no  matter  what  the  diagonal  elements  of  the  potential  matrix 

are  (attractive  or  repulsive),  if  the  off-diagonal  elements 

are  large  enough  in  magnitude,  there  will  be  a  bound  state^ 

because  only  the  square  of  g  enters  in  .  This 

corroborates  the  fact  that  the  coupling  term,  e.g.  Vl2_  ,  of 

9 

the  Hamiltonian  always  is  effectively  attractive 

However,  there  is  no  way  of  identifying  the  channel  in  which 

the  bound  state  occurs,  it  belongs  to  both  channels.  For 

a  bound  state  embedded  in  the  continuum  with  bound-state 

energy  E'  ,  0<CE!  <  80  ,  we  must  look  for  a  root  of  (3*23) 
b  b  c 

at  a  positive  real  value  of  k1(=  '/"eJ  )  and  at  a  positive 

b 

imaginary  value  of  k^(  —  4-i  ^E'  "  E2)J 
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p  2-ki -ie 
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cr  2 

g  1  2 


f  „ 

v2  (p  )p2dp  v2 (p ' )p  *  2 dp ' 


2-k?-ie 


J 


=  0 


p ' 2+k2 
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(3.26) 
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The  real  and  imaginary  parts  of  this  equation  must  vanish, 
i.e. 
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1  +  g2 


v2 (p )p  2dp 
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P2+k| 


=  0 
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Sl 1 [1+S22 
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v2  (p  )p  2dp-| 

P2+K2 


gl  2 
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v2 (p )p  2dp 


P2+K2 


J 


(3.27) 


The  first  equation  determines  E*  ,  so  g  should  be  negative 
(attractive),  while  the  second  equation  will  require  g^  =■  0 
and  g i l  may  be  positive  or  negative.  This  means  that 
channels  1  and  2  are  not  coupled,  then,  of  course, 
channel  2  can  have  bound  state  above  the  threshold  of  chan¬ 
nel  1.  However,  this  is  a  special  situation,  we  expect  that, 
for  more  than  two  channels,  some  zeros  of  A  would  give  rise 
to  bound  states  embedded  in  the  continuum  without  too  severe 
a  restriction  on  the  off-diagonal  elements  of  the  potential 
matrix. 


Resonant  poles,  in  the  one-channel  problem,  occur  in 


. 
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the  fourth  quadrant  close  to  the  positive  real  axis  in  the 

complex  momentum  plane.  In  the  two-channel  case,  there  are 

two  fourth  quadrants,  one  for  each  channel  momentum.  Poles, 

in  either  one  or  both  of  the  two  fourth  quadrants  and  close 

to  the  real  axis,  will  correspond  to  resonances.  When  the 

resonances  are  sharp,  i.e.  the  width  of  the  resonances  is 

small  compared  to  their  spacing,  a  more  convenient  way  of 

determining  the  resonance  energy  Er  is  to  require  the  real 

2.8 

part  of  A  to  vanish  at  E^  .  In  this  way,  we  are 

working  with  the  real  energy  E  instead  of  complex  E.  On 

writing  (3.23)  into  real  and  imaginary  parts, A(E)  =  A(e)+  iAT(E), 

R 

E^  is  then  determined  by  Ar  (E  )  =■  0  and  the  resonance 
width  F  is  proportional  to  Aj-(Ef).  For  E  near  Er  , 
we  can  expand 

A(E)  &  (E-E  )A'(E  )  +  iA  (E  )  (3-28) 

r  R  r  i  r  > 

where  A'r(e)  is  the  derivative  of  Ar(E)  with  respect 

to  E.  Thus,  near  resonance,  the  partial-wave  amplitude  a^ 
has  the  denominator  (3.28)  which  can  be  put  into  the  usual 
form  of  resonance  denominator. 

(E-E  )  +  iA  (E  ) / A ’ ( E  )  =  (E-E)  +  ir/2  ,  (3-29) 

r  i  r  r  r  j~ 


'  ' 

: 


where  P  is  the  width  of  the  resonance.  The  channel  in 
which  the  resonance  occurs  is  determined  by  the  conditions 
of  Ejp>0  or  .  If  Er<0,  then  Aj-=0  (since  A  is 

real  for  E<0)  so  that  P  =  0.  Hence,  there  is  no 
resonance  below  the  lowest  threshold  E  =  0.  We  can  rewrite 
a  (3.22)  with  the  common  resonance  denominator  (3.29)  for 


aag(E)~  ra6/[(E-Er)+ir/21 


(3-30) 


where 


r 


a£ 


(E  )/A'(E  )  . 
R 


(3-3D 


which  may  be  taken  to  be  real,  i.e.  principal  values  of  the 
integrals  in  n  .  are  taken,  since  the  width  V  is  assumed 

Ov  p> 

to  be  small.  If  0<Ej?<8^  ,  then,  from  (3  •  23 )  j  the  integrals 
involving  k,  (=  Er)  will  be  complex,  while  those  involving 

(=  Er  -  )  will  be  real.  We  can  find  F  from  Eqs.  (3.29), 

(3.23),  (3.22),  and  (3.31). 


r/2  =  A  (E  )/A'(E  )  =  -riiC  , 
I  R 


(3-32) 


a  A  30'  :)  0  s  t  » 0  >  •  1  •  c  1 

v 

( fS.£)  noil  ,n9ritf  ,  s3  >  3>0  'll 
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which  is  the  partial  width  for  channel  1.  If  E^>  £2  >  0, 
then  all  integrals  in  (3-23)  will  be  complex,  and  we  find 
similarly 

r/2  =  A  ( E  )  /A '  ( E  )  =  -rn/E"  -  r22/E  -e2  ,  (3-33) 

I  r  R 

which  shows  that  if  both  channels  are  open  (E  >£z),  then  the 
total  width  is  the  sum  of  the  partial  widths  for  each  channel. 


' 

, 
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Chapter  4.  A  SOLVABLE  MODEL  IN  NONRELATIVISTIC 

ELECTRODYNAMICS 

In  this  chapter  we  shall  discuss  the  interaction 
of  a  nonrelat ivistic  charged  particle  with  an  electromagnetic 
field.  The  charged  particle  is  assumed  to  be  spinless 
and  bound  in  a  nuclear  or  atomic  potentials  i.e.  it  satisfies 
the  ordinary  Schrodinger  equation.  An  example  is  furnished 
by  the  electron  of  a  hydrogen  atom  interacting  with  an 
electromegnet ic  field.  We  shall  limit  ourselves  to 
processes  involving  a  single  charged  particle  with  the 
emission  or  absorption  of  one  photon. 

After  a  description  of  the  model  (Section  1.), 
which  is  similar  to  the  well-known  Lee  model,  and  the 
discussion  of  the  processes  of  the  interaction,  we  arrive  at 
the  Hamiltonian  that  will  be  used  throughout  this  chapter. 

Our  main  interest  is  in  the  decay  behaviour  of  a  prepared 
state  (Section  2.)  which  leads  to  the  level  shift  and  the 
line  breadth  of  the  emitted  photon  (Section  3*)*  We  show  that 
the  prepared  state  is  physically  acceptable  (Section  4.). 

In  Section  5*  we  obtain  a  complete  set  of  states  and  give  an 
alternative  form  for  the  Hamiltonian,  since  the  results  of 
previous  sections  are  also  obtainable  by  expanding  the 


. 
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i 
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prepared  state  in  such  a  complete  set.  We  conclude  this 
chapter  by  considering  a  simple  form  of  spin-dependent 
interaction  in  the  Hamiltonian. 

1 .  The  Hamiltonian 


The  complete  Hamiltonian  for  a  nonrelativist ic 

bound  charged  particle  of  zero  spin  interacting  with  an  elec- 
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tromagnetic  field  is,  (fi  =  c  -  l) 


H 


com 


[(l/2m)Vi|j*(r,t)»Vi{»(r,t)+V(r)i|j*(r,t)iJj(r,t)]dr 
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(1/8tt  ) 
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+  (e2/2) 


r 


J ) 


^(rst)^?'  }t)ip(r’  ,t)ip(r,t) 


->  -> 
drdr  ’ 


r  -  r 


(4.1) 


The  first  term  describes  the  charged  particle  with  mass  m 

bound  in  a  nuclear  potential  V.  The  nucleus  can  be 

considered  fixed  (infinitely  massive).  The  second  term, 

where  we  have  used  Gaussian  units,  is  the  radiation 

energy  of  the  electromagnetic  field.  Both  E  and  A  are 

transverse  fields.  The  third  and  fourth  terms  give 

rise  to  the  interaction  between  the  charged  particle 

and  the  electromagnetic  field.  However,  the  fourth  term  is 

much  smaller  than  the  third  and  can  produce  two-photon 

processes.  The  last  term  is  the  electrostatic  interaction 

between  charges  which  does  not  contribute  in  one-electron 

processes.  Hence,  we  shall  omit  the  last  two  terms  in  the 
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Hamiltonian  for  our  discussion.  We  shall  come  back  to 

consider  the  fourth  term  in  the  next  chapter. 
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Let  us  expand  y  and  f *  in  a  complete  Prthonormal 

set  of  eigenfunctions  for  a  single  electron  in  a  nuclear 

-*■ 

potential.  The  vector  field  A  can  also  be  expanded  in 
plane  waves  with  periodic  boundary  conditions  in  a  box  of 
unit  normalization  volume.  We  then  have 


iKr,t)  =  Ibn(t)un(r) 
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where  u^r)  form  a  complete  orthonormal  set  of 
eigenfunctions  which  satisfy  the  Schrodinger  equation 

( -V2/2m  +  V(r))un(r)  =  Enun^  .  (4.3) 

The  photon  energy  is  =  (  k^l  and  e^  is  the  unit 

polarization  vector  orthogonal  to  k^  .  The  sum  over  A_ 
includes  momentum  eigenstates  and  polarization,,  different 
direction  of  polarization  is  denoted  by  different  k 
The  expansion  coefficient  b^  (bn  )  is  the  creation 
(annihilation)  operator  for  an  electron  in  the  electron 
eigenstate  n.  Similarly,  a  A.  (aA.  ^  creates  (annihilates) 
a  photon  of  momentum  ,  polarization  .  The  b's  and  a's 

satisfy  the  usual  commutation  relations. 
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[ax(t),ax,(t)J  =  [aj(t),aj.,(t)]  =  0 


[aA(t)’bn(t)1  =  [ax(t),b+(t)] 


[aj(t),bn(t)]  =[a+(t),b+(t)] 


(4.4) 


We  farther  assume  that  there  are  only  two  energy  levels  in 
(^•3),  or  only  two  eigenstates  contribute  significantly  to 
the  interaction,,  the  higher  level  Eg  and  the  lower  level 
Ej  ,  e.g.  we  can  take  two  bound  states  in  the  hydrogen  atom. 
If  we  are  restricted  to  processes  of  the  following  type: 

One  electron  in  level  2*->one  electron  in  level  1  ■+-  one  photon 

(4.5) 

as  indicated  at  the  beginning  of  this  chapter,  we  will 
obtain  a  simple  Hamiltonian  by  substituting  (4.2)  into  (4.1), 
(ignoring  last  two  terms  in  (4.1)  ) 
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where 


t 
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H0  =  E2b2(t)b2(t)  +  Eibi(t)bi(t)  + 


Ik 

A 


xaJ(t )aA(t ) 


and 

Hi  =  e^(l//2k^)(0xb|(t)bi(t)ax(t)+3«bT(t)aJ(t)b2(t))  .  (4.6) 

The  infinite  zero-point  energy  of  the  radiation  field,  2  k^  /2? 
has  been  subtracted.  The  quantities  (3^  are  given  by 


3  x  =  ( i  /47r/m ) 
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ik  *r 

uf(r)e  ?L 
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Vui (r )dr 


(4.7) 


j 


which  are  essentially  the  one-photon  interaction  matrix 
elements  (p.l43,  Heitler)^1. 

The  Hamiltonian,  as  given  in  (4.6),  is  in  Heisenberg 
representation  since  the  a's  and  b's  are  time-dependent 
operators.  We  can  transform  (4.6)  into  Schrodinger 
representation  which  coincides  with  Heisenberg 
representation  at  a  particular  time  t  =  0.  The  transforma¬ 
tion  of  operators  is  achieved  by 
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0(t)  =  elHt  O(t=0 )  e-iHt  . 


(4.8) 


Hence,  the  Hamiltonian  in  Schrodinger  representation  has 
the  same  form  as  (4.6)  but  with  the  a's  and  b's  specified 
at  t  =  0  which  we  shall  simple  rewrite  omitting  the  time 
parameter.  In  the  following  sections  the  Hamiltonian  in 
Schrodinger  representation  will  be  used,  since  we  shall 
try  to  solve  the  time-dependent  Schrodinger  equation  for  a 
prepared  state  of  our  physical  system.  The  Hamiltonian 
(4.6)  is  very  similar  to  that  of  the  Lee  model  5  , 

except  that  we  do  not  have  to  introduce  a  "cutoff"  for  the 
interaction  energy,  since  ,  as  defined  in  Eq .  (4.7)  will 
provide  a  natural  cutoff  for  all  integrals  involved. 

So  far  we  have  not  really  specified  the  two  levels 
of  the  bound  electron  (4.3).  If  those  two  levels  are  two 
discrete  states  of  an  electron  in  an  atom,  then  Hamiltonian 
(4.6)  describes  radiative  transitions  between  two  discrete 
states  by  emission  or  absorption  of  one  photon.  Generali¬ 
zation  to  transitions  of  an  atom,  or  nucleus,  from  any 
state  of  higher  to  one  of  lower  energy  can  be  easily  carried 
out.  We  can  choose  the  appropriate  wave  functions 
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for  those  two  states,,  either  or  both  of  which  may  be  in  the 
discrete  or  in  the  continuous  energy  spectrum,,  Hence,  for 
atoms,  the  photoeffect,  the  radiative  capture,  (discrete 
continuous  transitions)  and  Bremsstrahlung  (continuous 
continuous  transitions)  can  be  studied  along  the  same 
lines.  Similarly  for  nuclear  photoeffect  and  Y_(^ecay 
of  nuclear  levels,  but  the  wave  functions  for  the  charged 
particle  in  a  nuclear  field  (not  necessarily  Coulomb 
potential)  will  be  very  complicated  or  simply  unobtainable 
since  the  behaviour  of  nuclear  forces  for  very  short  ranges 
is  not  known. 

2 .  The  Development  of  a  Prepared  State 

We  will  be  using  Hamiltonian  (4.6)  in  Schrodinger 
representation.  Let  us  first  define  a  vacuum  state  |0^> 
for  our  system,  i.e.  a  state  with  no  electron  and  photon 
present. 


b  2  |  0  >  =  b i  | 0  >  =  a^|0>  =  0  ,  (4.9) 

which  is  also  an  eigenstate  of  (4.6)  with  eigenvalue  zero. 
This  vacuum  state  is,  therefore,  time-independent  and  we 
can  build  up  the  states  we  want  for  our  system  by  applying 
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a 1 s  and  b '  s  to  it . 

Suppose  that  we  prepare  our  system  at  t  =■  0 
to  be  at  the  higher  level  2  with  no  photon,  i.e.  in  the 
initial  state 

¥(0)  =  bj  |0>  ,  (4,10) 

and  we  shall  study  its  subsequent  development  due  to  the 
Hamiltonian  (4.6).  At  a  later  time  t  >  0,  our  system 

will  be  in  the  state  \{/(t)  which  is  governed  by  the  time- 
dependent  Schrodinger  equation, 

i3H'(t)/at  =  H'F(t)  ,  (4.11) 

where  H  is  the  Hamiltonian  operator  (4.6).  The  formal 
solution  of  the  time-dependent  wave  function  ^(t)  has  the 
explicit  form 
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(4.12) 
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which  satisfies  the  specified  boundary  condition  (4.10) 

42 

at  t  =  0.  For  t  ?  0,  we  can  rewrite  (4.12)  as 
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¥(t)  =  (  ~  1/  2  7T  i  ) 


dEe"iEt 

FT? - u —  Y(0), 

E+ie-H  ’ 


loo 


(4.13) 


where  £-^  0  (+)  and  E  is  a  real  variable. 


Our  main  problem  now  is  to  find  the  quantity 


E+ie-H  ^  E+ie-H  b2l0>  * 


(4.14) 


Using  the  identity,  valid  for  operators  as  well  as  c-numbers. 
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and  putting  a  =  E  +  i£  -  H0  ,  b  ^  -H1  ,  we  can  write 

(4.14)  as 


1 

E+ie-H 


1 

E+ie-E2 


0>  +  el 

X 


(4.16) 


Using  a  formula  in  Schweber  p. 
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E+ie-H  E+ie-H-k.  ’ 


(4.17) 


we  obtain 
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Substituting  (4.18)  into  (4.16),  we  finally  get 
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(4.19) 


Eq.  (4.13)  together  with  the  above  expression  determines 
the  complete  behaviour  of  our  system  which  was  prepared  in 
the  state  "''p'(O)  at  t  =  0. 

We  are  mostly  interested  in  the  probability  of 
finding  the  system  in  its  initial  state  after  time  t.  The 
probability  amplitude  for  the  system  remaining  in  ^(O) 
at  a  later  time  t  is  given  by 
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(4.20) 


In  the  denominator.,  the  sum  over  \  (including  sum  over  two 
directions  of  polarization)  can  be  written  as  an  integral 
when  the  normalization  volume  tends  to  infinity. 
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which  is  defined  for  real  E.  We  have  assumed  that  does 

— ^ 

not  depend  on  the  direction  of  the  momentum  vector  but 

only  on  its  magnitude  Ik^l  .  The  symbol  P  stands  for 
the  principal  value  of  the  integral.  Let  us  write  down  the 
denominator  in  (4.20)  as 

hCE-Ei-ie)  =  E+ie-E2+f ( E-Ei +ie )  5  (4.22) 


which  is  also  defined  for  real  E,  and  can  be  analytically 
continued  for  a  general  complex  E  provided  that  the  same  can 

2_  3  9 

be  done  for  |(3>(k)| 

Let  us  define  a  complex  variable 


z=E-Ei=x+iy  j  (4.23) 

where  E  is  now  complex,  and  E4  is  real,  the  lower  energy 
level.  On  the  first  Riemann  sheet,  for  0  <  arg  (z) <  2 7T 


1 


h~^(z)  =  z  +  Ei  -  E2  +  f(z) 


(4.24) 
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which  has  a  cut  along  the  positive  real  axis  as  can  be  easily 
seen  from  (4.21)  and  (4.22).  Integral  (4.20)  may  be  obtained 
by  the  following  contour  integration  on  the  first  Riemann  sheet: 


Fig.  1. 


Contour  for  Integral  (4.20)  Closed  on  the  First  Sheet. 


' 


We,  therefore,  must  find  the  roots  of  h1  (z)  in  the  lower 
half  plane,  excluding  the  cut.  A  simple  substitution  of 
z  —  x  4-  iy  into  h^  (z)  =  0  will  require 
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y  =  0 


and 
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x  =  E2-Ei-(e2/27T2 ) 


3(k)  I 2kdk/(k-x) 


(4.25) 
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for  x  <  0.  Eq .  (4.25)  has  no  solution  below  x  =  0,  if 


(4.26) 


which  we  shall  assume  and  E±  satisfy.  We  expect  that 

the  behaviour  of  (4.20)  is  some  kind  of  damped  exponential 
which  arises  from  a  complex  root  of  the  denominator  with 
negative  imaginary  part.  Hence,  we  search  further  onto  the 
second  Riemann  sheet. 

Analytic  continuation  onto  the  second  Riemann  sheet, 
for  0  >  arg  (z)  >  -2 71  ,  is  effected  by  defining  for 

real  z  (  =  x) , 


h11 (x-ie)  =  hI(x+ie)  . 


(4.27) 
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Then,  for  z  x  +  iy,  we  have^ 

h^(z)  =  h^(z)  +  i  ( e  2  /tt  )  |  3  ( z  )  |  2  z  . 

Now  we  can  evaluate  Integral  (4.20)  by  closing  the 
contour  on  the  second  Riemann  sheet  as  shown  in  Fig. 


z  plane 
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(4.28) 
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2. 


Fig.  2. 


Contour  for  Integral  (4.20)  Closed  on  the  Second  Sheet. 
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The  solid  portion  of  the  contour  is  on  the  first  Riemann 
sheet  while  the  dotted  portion  is  on  the  second.  The 
integral  of  the  closed  contour  can  be  obtained  from 

Cauchy's  integral  formula,  so  we  must  find  the  zeros  in 
h11  (z).  Assuming  for  simplicity  that  h11  (z)  has  only  one 

solution,  i.e.  the  integrand  in  (4.20)  has  a  simple  pole,  at 

z  -  (E-^  ~E  1  )  -  i  Y/2,  then  from  (4.28),  (4.24),  and 

, .  .  38 

(4.21)  we  obtain  the  determining  equations  for  E  and  T, 
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where 


C  +  in  =  |3(Er-E1-iy/2) |2(Er-E1-iy/2)  .  (4.30) 


Integral  (4.20)  can  be  written  as  contour  integrals 


<Y(0)  |'F(t)> 
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where  C1  and  C2  are  shown  in  Fig.  2.  The  denominator 
h  (z)  will  be  either  hr  (z)  or  h11  (z)  depending  on 
whether  the  contour  is  on  the  first  or  on  the  second 
Riemann  sheet.  The  integrals  over  C±  and  are  then. 
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and 
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(4.32) 


where  we  have  used  (4.28)  with  z  =.  x  (real)  for  the 
integral  over  C2  .  On  comparing  the  above  two  integrals, 
we  see  that  for  weak  coupling  and  not  too  large  t  the  first 
one  dominates  while  the  second  one  will  only  become 
appreciable  after  a  sufficiently  long  time.  However,  by 
this  time  (4.3l)  may  have  already  become  extremely 
small  indeed.  As  t  tends  to  infinity  the  behaviour  of  the 
second  integral  in  (4.32)  is  determined  by  the  "critical 
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point  at  x  =  0  ,  which  means  that  most  contribution 

comes  from  the  end  point  x  =  0.  We  can  obtain  the 
asymptotic  expansion  for  the  second  integral  by  expanding 
the  square  bracket  about  x  =.  0, 
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B(x) 


hI(x)hII(x) 
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]  n+  x^— [ 
x=0  dx 


1 x=0  +  ' 


(4.33) 


and  by  inserting  the  above  expansion  into  the  integral 
which  can  then  be  evaluated  at  x  =  0  to  yield 
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The  manner  in  which  (4.34)  is  obtained  tells 

us  that  the  asymptotic  behaviour  of  the  probalility 

amplitude  depends  on  the  dynamic  structure  of  the  system 

40 

and  on  how  one  prepares  the  initial  state  .  For  example,, 
if  we  start  with  the  system  in  the  initial  state 


Y(0) 


(4.35) 


which  describes  an  electron  in  level  1  plus  a  packet 
of  photons  of  arbitrary  shape  .  A  similar  calculation 

will  give  the  exponential  form  for  small  t,  and  for  large  t 
an  integral  over  „ which  is  dominated  by  the  contribution 

from  x  =  0, 
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Again,  we  expand  the  square  bracket  about  x  =  0,  like 

(4.33) 
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where  tj>^_  is  the  arbitrary  shape  of  the  wave  packet  that  we 
prepare  our  system  at  t  =  0.  This  means  that  we  can  pick 
the  leading  term  in  t  for  the  asymptotic  expansion  of 
(4.36)  by  choosing  suitable  <j>^  ;  e.g.  if  the  first  term 

in  (4.37)  vanishes  for  a  particular  choice  of  cf>^  , 
then 


Eq . 


(4.36) 


constant(e  iEltt  +  0(t 


(4.38) 


-3 


-2 


which  has  a  leading  term  t 
(4.34) . 
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Let  us  come  back  to  the  time  evolution  of  our 
system.  The  probability  of  finding  the  system  in  the 

z 

initial  state  at  time  t  is  represented  by  |<  "ty(  o)  |  'ty(t)  I 
We  can,  therefore,  conclude  from  (4.3l).  (4.32).  and 
(4.34)  that  this  probability  decreases  exponentially  with 
lifetime  1/ Y  f°r  small  t  ( 7^  1/ Y  ),  and  ultimately 
decays  as  t”4  .  However,  one  must  note  that  the  preparation 

of  the  initial  state  affects  greatly  the  asymptotic  behaviour 
of  the  above  mentioned  probability. 

3.  Level  Shift  and  Line  Breadth 


The  complex  root  z  =  -  E L  -  i Y/z  for  hrI  (z) 

as  determined  by  (4.29)  has  simple  physical  interpretation. 

2. 

Let  us  first  look  at  the  weak  coupling  limit  in  which  e 
is  very  small,  E  r  &  E  z  ,  and  E  r  -  E  ^  Y/Z  2^.  0( ^ )  • 
In  this  limit  we  can  solve  for  Ep  and  T  without  knowing 
the  exact  form  of  (k) .  From  (4.29)  and  (4.30)  then. 
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and  therefore, 


3 ( k) | 2kdk 

k+E i — E  2 
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Y  »  (e2/ir)  |6(E2-E1)  |2(E2-Ei)  . 


(A. 39) 


The  results  for  dE  and  y  agree  with  that  by  first  order 

time-dependent  perturbation  calculation  for  the  level  shift 

11( §18) 

and  the  line  breadth  respectively.  When  the  electron 

at  level  2  decays  spontaneously  to  level  1  with  the 
emission  of  a  photon,  the  spectral  line  of  the  emitted 
photon  will  have  a  maximum  at  frequency  =  E  ^  -  Ej,  with  a 
half  width  Y  ,  instead  of  a  perfectly  sharp  line  with 
frequency^E^  -  E  j  .  Or  one  may  say  that  when  the 
electron  at  level  2  interacts  with  an  electromagnetic  field 
according  to  (4.6),  its  original  level  E^  will  shift  to  E^ 
with  a  breadth  Y  ,  i.e.  it  has  lifetime  of  l/y  due  to  the 
radiative  transition  to  the  lower  level  1.  Eq .  (4.39)  also 
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gives  the  correct  behaviour  of  our  system  when  there 

is  no  interaction  (e2  -  0),  such  that  aE  =  Y  =  0.  In 

this  case  the  electron  will  remain  at  level  2  indefinitely, 

„  2_ 
since  it  is  an  eigenstate  of  the  Hamiltonian  when  e  =•  0. 

Suppose  we  can  increase  the  coupling  strength 

2. 

e  between  the  bound  electron  and  the  electromagnetic 
field  without  affecting  the  other  part  of  our  system, 
i.e.  E  £  ,  E ^  ,  and  3(k)  will  not  change  if  e2-  changes. 

A  look  at  the  second  equation  of  (4.29)  will  show  us  that 
there  is  another  solution  for  with  e  ^  ^  0. 

o  2. 

If  such  a  solution  exists  for  a  certain  e  -  e s  , 
then  we  would  have  an  emitted  photon  with  a  perfectly  sharp 
spectral  line.  This  solution  will  require  £  -  0  which 
implies  Er  =  E  ^  When  this  condition  is  applied  to 
the  first  equation  of  (4.29),  we  find  that  E^  and  E^ 
must  satisfy 
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which  is  the  equation  that  determines  e2  .  However,  in 
order  to  obtain  (4.29)  such  that  the  electron  at  level  2 
can  decay  spontaneously,  we  had  already  assumed  that  E1  and 
E  ^  must  fulfil  the  inequality  (4.26)  which  holds  for  all 
values  of  e'1  .  Eq .  (4.40)  is,  therefore,  invalid,  and 

there  is  no  emitted  photon  with  perfectly  sharp  line.  Let 
us  suppose  for  the  moment  that  e2-  e^  without  the  above 

contradiction,  can  e2  be  increased  further  still?  If 
Er  and  T  are  both  continuous  functions  of  ez  ,  when  e2 
increases  from  0  to  es  ,  E^,  decreases  from  E2  to  Ex 

and  T  from  0  back  to  0.  So  if  e2  >  ej  ,  E^  =  E1  -  J  E 
where  <?E  is  a  small  positive  energy.  Substitution  of 
Ep  -  E^  J'E  into  the  second  equation  of  (4.29)  will 
yield  a  negative  Y  which  means  <  ¥(o)  |  ^(t)}  grows  indefinite¬ 
ly,  and  hence  not  acceptable,, 

We  conclude  this  section  by  remarking  that  we 

O 

cannot  increase  the  coupling  parameter  e  infinitely  even 
if  it  does  not  affect  the  other  structure  of  the  system. 

The  exact  form  of  §  (k)  must  be  specified  if  we  want  to 
find  the  exact  solutions  for  Er  and  Y  . 
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The  Mean  Energy  and  the  Rate  of  Energy  Dissipation  for 

the  Initial  State 


We  have  prepared  our  system  at  t  =■  0  in  the 
initial  state  (4.10)  which  evolves  according  to  the  time- 
dependent  Schrodinger  equation  (4.1l).  We  have  also 
studied  the  decay  behaviour  of  the  system  in  Section  2. 

In  this  section,  we  want  to  show  that  the  initial  state  (4.10) 
for  our  system  is  physically  acceptable. 

The  mean  energy  and  the  rate  of  energy  dissipation 
for  an  acceptable  state  should  be  expected  to  be  finite. 
Otherwise  we  cannot  simply  prepare  at  t  =  0  such  a  state 
of  infinite  energy  or  a  state  of  finite  energy  but  with 
infinite  time  derivative.  In  the  latter  case  the  initial 
state  jumps  discontinuously  into  another  state  and  thus 

(t)  is  no  longer  a  continuous  function  of  t  as  described 
by  (4.1l).  The  mean  energy  of  the  initial  state  Sjf(o)  is 
determined  by  <  *<o)|HIY(0)>  where  \J/  (0)  and  H  are 
given  by  (4.10)  and  (4.6)  respectively.  A  simple  calculation 
will  yield  the  result 


<¥(0) I H | ¥ ( 0 ) >  =  E2  , 


(4.41) 


■ 


>f  r 


_ 
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which  is  just  the  bound  state  energy  of  the  electron  at  level 
2  and  has  a  finite  value.  The  rate  of  energy  dissipation  for 
2^  (0)  at  t  =  0  is  proportional  to  the  expectation  value 
of  in  the  initial  state. 


<¥(0)  |H2  |V(0)>  =  (E2)2  +  e2£|f3x|2/(2kA)  ,  (4.42) 

A 

by  a  similar  calculation  as  (4.4l).  The  first  term  in  (4.42) 
is  again  finite.  The  second  term  is  also  finite  with  (3^_ 
given  by  (4.7)  which  will  provide  the  necessary  converging 
factor  in  the  sum.  We  have  demonstrated  that  ^(o)  for  our 
system  is  physically  acceptable  and  its  subsequent  decay  is 
exponential  for  small  t,  and  for  large  t  asymptotically 
in  inverse  powers  of  t. 

We  should  mention  here  that  all  results  obtained 
so  far,  (4.20)  in  particular,  can  also  be  derived  by 
expanding  the  initial  state  in  terms  of  a  certain  complete 
set  of  eigenstates  of  the  Hamiltonian  (4.6),  providing 

and  E  ^  satisfy  condition  (4.26).  In  the  next  section 
we  shall  obtain  this  complete  set  of  eigenstates. 


suIkv  ioidfsJo9c  C9  'rii  ('6^  Isnc-i  jf^ocjo^q  Bx  0  =  3 


£ 

■ 


■ 
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5.  A  Complete  Orthonormal  Set  of  Eigenstates  of  the 

Hamiltonian  and  An  Alternative  Form  for  the  Hamiltonian 


We  construct  here  a  certain  complete  orthonormal 

set  of  eigenstates  of  the  Hamiltonian  (4.6)  in  the  subspace 

i  i  i  37 

spanned  by  b^|Q>  and  bj  a^|0>  .  These  eigenstates 

are  solutions  of  the  Schrodinger  equation 

H 1 1 ,  A>  =  (Ei+kx)|l5A>  ,  (4.43) 

where  II,  A.')  are  the  eigenstates  which  are  states  describing 
the  scattering  of  a  photon  by  the  electron  in  level  1.  we 
can  write  as  the  electron  wave  function  in  level  1 

plus  an  incoming  plane  wave  of  a  photon,  bj  a^|0>  ,  and  an 

outgoing  scattered  wave  ix>  ,  then 

| 1 , A>  =  btax|0>  +  | x>  •  (4.44) 

Substitution  of  (4.6)  and  (4.44)  into  (4.43)  yields  the 
solution  for  the  outgoing  scattered  wave 


egA 

/2iy 


Ei +k, -H+ie 


-hi 


A 


X> 


0> 


(4.45) 


' 
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. 


70 


-1  X 

The  operation  of  (Ej_  +  -  H+  i£  )  on  b^|o> 

is  given  by  (4.19)*  so  we  obtain  finally. 


|l,A>  =  btat |0> 


e  3 


A 


A 


/2k 


A 


(Ei+k,+ie-E2+e2y 
A  A' 


Is 


1 


-1 


2k 


A 


kA'“kA"le 


( b  2  |  0> 


6*„  btaj„|0> 

el  —=  - )  • 

A  M/2k^„k^  ,,-k^-ie 


(4.46) 


The  proof  that  form  a  complete  orthonormal  set  in  the 

subspace  mentioned  at  the  beginning,  namely 


J  1 1 ,  A><A  ,1 1  =  Jbta|  |  0x0  |  a,  b  !  +  bj|0><0|b2 
A  A 


<A,1|1,X'  >  =  6aa, 


(4.47) 


' 
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37 

is  the  same  as  given  by  Glaser  and  Kallen  ,  where  we 
have  the  occasion  to  evoke  the  inequality  (4.26). 

Next  we  shall  write  an  alternative  form  for  the 

Hamiltonian  which  is  quadratic  in  oprators  and  is  very 

45  J 

similar  to  Henley’s  .  We  note  in  the  subspace  that  we 

have  considered  the  Hamiltonian  will  be  operating  on  bj  |  o> 
and  bJl(X>  (plus  a  packet  of  photon).  From  the  type  of 
allowed  transitions  (4.5)*  we  would  like  to  have  operators 
that  change  bj,jO)>  to  bJiOy'  and  vice  versa.  Let  us  define 
new  operators  T+  and  T_  =  T  +  that  possess  the  desired 
property, 

T+  =  b|b2  ,  and  T_  =  bjbx  ,  (4.48) 


such  that 


t 

T+Cb  2  |  0  >  ) 


T+(bt|0>)  =  0  , 


T  (bt |0>)  =  0 


T 


(bt 


0>)  = 


bt 


0> 


(4. 49) 


. 
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The  commutator  of  T+  and  T_  is  defined  as 


T3  —  [ T+ ,  T_ ]  —  bib 1  —  b|b 2  » 


(4.50) 


which  has  the  following  property 


T  3 ( b  2  I  0>  )  =  -bj  I  0  >  3  T  3 ( b t I 0> )  =  bt  I  0> 


(4.51) 


In  our  particular  subspace  spanned  by  b!>|0)> 
the  operator  (b|  b^  +  b£  b^  )  has  eigenvalue 
write 


and  bj  aj^|o^> 
1,  so  we  can 


bjb2  =  (l-Ts)/2  ,  btb,  =  (l+Ts)/2  .  (4.52) 


Substitution  of  (4.48)  and  (4.52)  into  (4.6)  yields  the 
result  for  the  Hamiltonian 


H  =  (Ez+EO/2 


(E2-E1)T3/2  + 


A 


AaAaA 


+  eZ(2kx)Js  (gxaxT_  +  3*a+T+) 

A 


(4.53) 


which  is  meaningful  only  in  the  special  subspace  concerned. 


' 

mmi  «  -  dU  ■  i*  ^  =  ’  % 

; 

.  ,  4  .  <0  j  d  :<0J'  <t*|  t<5-  =  (<0U<J)tT  *' 

■  '  .  .•  '  '  '  '  ' 

.  (et+-:  )  i  -  •  5  t  :  \  ( ,  = 
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6 .  Spin-Dependent  Interaction 

The  interaction  between  electron  spin  and  the 
radiation  field  will  modify  the  Hamiltonian  (4.1),  and  an 
additional  term. 


H 

sp 


( -e/2m) {  i|;*(r,t)a‘H(r,t)i|4r,t)dr 


=  (-e/2m)  1^* (r ,t )a • ( V*A(r ,t ) ) ty(r ,t )dr 


(4.54) 


should  be  added  to  (4,l).  In  this  equation  cr  is  the  usual 
Pauli  spin  matrix,.  The  spin-orbit  interaction  is  neglected 
in  the  nonrelativistic  limit.  However,  the  contribution  from 

H^r.  is  still  extremely  small  in  causing  transitions  when 

r 

compared  with  the  third  term  in  (4.l),  unless  the  third  term 
is  highly  forbidden  for  a  certain  kind  of  transition.  We  are 
interested  in  HCK  mainly  because  it  is  linear  in  A  and 
thus  involves  processes  of  one  photon  which  may  be  combined 
with  the  third  term. 


A  similar  expansion  for  the  first  three  terms  in 


br 


" 
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(4.1)  and  can  be  carried  out  as  before  but  with  V/ 

now  being  a  two-component  field. 


p(r,t ) 


I 


/b n+<*>\ 

b  ( t )  ) 

n-  / 


u  (r ) 
n 


(4.55) 


where  -  u^(r)  are  eigenfunctions  of  (4.3).  Operators  bm 
and  bn_  are  annihilation  operators  for  an  electron  in  state 
n  with  spin  up  and  spin  down  respectively.  For  simplicity, 
the  magnetic  field  H  is  assumed  to  be  polarized  in  the 

— y 

z  direction,  i.e.  the  electromagnetic  vector  potential  A 
is  propagated  in  the  x  direction  and  polarized  in  the  y 
direction,  so  we  obtain 


H(r  ,t ) 


ik 


=  V*A(r,t)  =  £i/2iTk^z(a^  (t  )e 


X 


-> 
'  V 


-aj^  ( t )  e 


-lie 


A 


■r 


(4.56) 


where  z  is  a  unit  vector  in  the  z  direction.  Note  that 
sum  over  X  here,  and  what  follows  in  this  section,  sums  over 
momentum  eigenstates  only,  since  there  is  only  one  direction 
of  polarization,  the  y  direction.  Then 


' 

■ 
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?• ( VxA) 


Ji/27rk^  (a^(t  )e 
A 


->  ->  ->  -> 
ik, • r  ,  -ik,  *r 

A  -  a+(t)e  A  ) 


1  0 

0  -1 


(4.57) 


If  we  again  consider  a  two-level  system  with  allowed  transi¬ 
tions  of  the  type  (4.5)*  a  similar  Hamiltonian  as  (4.6)  can 
be  obtained  in  Schrddinger  representation, 


H  =  E2(bJ+b2+  +  b|_b2_)  + 


Ei(bt+bi+  +  b t_b i _ )  + 


Jk,  ata. 


A 


"AAA 


+  e£  ( 1/ /2k ^ ) (3^(b2+bi+ 
A 


t 

b  I 


-bl-)aA  +  B*(bl+b2+  +  bt 


b  ■ 


)a+) 


+  e  J  ( 1/  /2k^)  (  ot^  (  b  2+b  i  + 


a*(bl+b2+ 


bt 


b 


)  aA ) 


(4.58) 
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where  (3-^  are  given  by  (4.7)*  and 


lit 


aX  = 


( -i/frk^/m) 


luf  (r  )e 


X 


-> 
'  T 


*  y 

ui (r)dr 


(4.59) 


Operators  a^  and  bn+.  obey  the  commutation  relations  (4.4) 
with  the  addition  that  operators  of  different  spin  commute. 

4- 

The  development  of  a  prepared  state  like  (4.10),  M7(o) =  I oy  , 
can  be  analysed  following  the  same  procedures  as  before  and 
similar  results  obtained.  For  example,  the  probability 

amplitude  at  time  t  of  finding  the  system  in  the  initial 

1 

state  b£+  |  0^>  is  given  by 


oO 

r 


<T(0) |¥(t)> 


-1 

2-ni 


dEe 


-iEt 


J 


E+ic-E2  +e  2  J- 
X 


3X+aX 


1 


2k 


X 


k. +Ei -E-ie 
X 


(4.60) 


which  can  be  compared  with  (4.20)  and  discussed  in  the  same 


way . 
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Chapter  5.  APPROXIMATION  FOR  THE  Aa  TERM  IN  THE 

HAMILTONIAN 

It  was  mentioned  in  Section  1  of  last  chapter 
that  we  shall  return  to  study  the  Az  term  in  the  Hamiltonian 
(4.1).  We  can  see  from  expansion  (4.2)  for  kz  that  this 
term  will  connect  states,  initial,  final,  or  intermediate, 
that  differ  by  two  or  no  photons.  In  Section  1,  we  make  the 
approximation,  valla  for  small  t,  which  is  equivalent  to 
fixing  the  number  of  electrons  at  t  =  0  In  our  system. 

This  approximate  Ha  Is  then  treated  together  with  Hraj  , 
Eq.  ( 5 • 7 ) j  classically  by  canonical  transformations  to  yield 
a  diagonalized  form  which,  in  turn,  is  quantized  (Appendix). 
The  quantized  expression  (A. 13)  contains  a  new  zero-point 
energy  which  Is  different  from  the  original  zero-point 
energy  of  the  radiation  field.  Only  when  the  coupling 
strength  ez  vanishes  identically,  are  these  two  zero-point 
energies  equal,  since  H^  will  also  vanish  identically. 

The  old  creation  or  annihilation  operator  for  the  photon  will 
be  transformed  into  an  infinite  sum  of  new  creation  and 
annihilation  operators.  This  implies  that  a  "free"  photon 
corresponds  to  a.  super^oj/h^K?  f  "  transformed  "  photons  (Section  2). 

Section  3  deals  with  the  “physical”  and  the  "bare"  vacua 
for  the  photon  and  the  possibility  of  preparing  our  system 
at  t  -  0  in  terms  of  the  bare  vacuum.  If  one  can  switch 


. 
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off  initially,,  one  is  led  strongly  to  the  idea  that  an 

initial  state  can  be  built  up  from  the  bare  vacuum.  Under 
very  special  assumptions,  that  both  (5. 15)  and  (5. 19)  are  fi¬ 
nite,  this  idea  is  indeed  possible.  However,  we  would  rather 
take  the  usual  notion  in  quantum  field  theory  that  one  should 
not  and.  cannot  separate  out  a  bare  particle  from  the  physical 
particles.  So  a  physically  acceptable  state  can  only  be 
prepared  from  the  physical  vacuum. 


The  Az  term  in  (4.l)  is  given  by 


(e2/2m) 


A2(r,t)b*(r,t)b(r,t) dr  , 


(5.D 


where  A  is  the  electromagnetic  field  and  Y'  the  field  for 
the  bound  electron.  Our  approximation  is  to  replace 


|iKr,t)|2^>  k(?,0)|2  ,  (5.2) 

which  is  no  longer  the  dynamical  operator  but  a  real  form 
factor,  e.g.  the  charge  distribution  for  the  bound  electron 
at  t  —  0.  This  approximation  is  valid  for  small  t  and  we 
shall  be  mainly  interested  in  studying  the  behaviour  of  our 
system  at  small  t.  The  charge  distribution  may  be  normalized. 
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i .  e  „ 


r 

|iK?,0) | 2dr  =  1  ,  (5-3) 

if  we  assume  that  at  t  -  0  there  is  one  electron  in  our 
system,  at  level  2  or  1. 

For  the  time  being,  we  shall  treat  A  as  a  classical 
field  and  make  the  following  plane-wave  expansion: 


J(r,t)  =  2/fFj6^(q^(t)cos(S^*r)  -  ( 1/k^  )p^  ( t )  sin(  k^  *r )  )  . 

X 

(5.4) 

The  amplitudes  q.  and  p  are  pairs  of  canonical  variables 
which  are  real0  Quantization  of  field  A  can  be  accomplished 

by  interpreting  p  and  q  as  operators  in  Heisenberg 

A-  x_ 

representation  satisfying  the  usual  commutation  relation. 


[qx(t) ,  PA , ( t ) 3  =  i«AA , 


(5.5) 


Further,  a  system  of  creation  and  annihilation  operators, 
a£  and  a^  ,  can  be  introduced  (second  quantization). 


LX 


=  /l/2k,  (a,  + 


X 


'X 


ts 

a,  ) 


P 


X 


=  i/kx/2 


( aA  -  aA> 


(5.6) 
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The  above  quantization  procedures  have  already  been  carried 

—  ^ 

out  in  the  previous  plane-wave  expansion  for  A  in  (4.2). 
We  use  expansion  (5.4)  here  because  we  want  to  treat 


f 


H 


rad 


+  H2  =  ( 1/8-rr ) 


->• ,  -> 


-> 


(E2(r,t)+(VxA(r,t) )2)dr 


+  (e2/2m) 


(r,t) |^(?,0) 


;dr 


(5.7) 


as  a  classical  Hamiltonian  expressed  in  canonical  variables 
p^  and  q^  .  After  H^^  +  H^  has  been  diagonalized  by 
canonical  transformations,  which  will  then  be  quantized  as 
indicated  in  (5*5)  and  (5.6),  introducing  a  new  set  of 
creation  and  annihilation  operators.  The  diagonalization  of 
(5.7)  is  shown  in  the  appendix. 

2.  The  Transformed  Hamiltonian 


In  this  "transformed"  Hamiltonian,  we  must  put 
the  previously  subtracted  infinite  zero-point  energy  of  the 
radiation  field,  ^  k^/2,  back  to  (4.6),  since  it  is  included 

r-\  v-'i  s-\  s~\  -P  /~s  /''v 

rad 


in  H„,j  ,  and  of  course  Hg  will  also  be  present.  The 


Hamiltonian  becomes 


b£91 

tfns 
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H 


T 


(H+Jk  /2)  +  H2 
A  A 


—  (E2b2b2  +  E i b i b i  +  H 


rad 


+  Ei)  +  H; 


(5.8) 


where  H  is  given  by  (4.6).  Using  Eqs.  (5.6),  (A. 3),  and 
(A. 8),  the  old  annihilation  (or  creation)  operator  a^  is 
transformed  into. 


a,  =  /kX72  I  /T72ir(a;(lvA-^vA)+a’  +  (yvX+yvA))  . 

V 

(5.9) 


where 


^  vA  ^/^AA,  +  (l/kA)^AA’  ^  (T  }vA' 


and 


^vA  ^/lKv^AA,+^l/kA)^AA’ )TAtv 

A 


(5.10) 
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The  new  set  of  annihilation  and  creation  operators,  a^ 
j  are  defined  by  relations  similar  to  (5.6).  With 
we  can  write 


HT  =  HT-hv/2  = 

V 


+ 


+1 


E2b2b2+  E^bjbjHi"^ 


v 


K  a,+ 
v  v 


a 


t 

v 


+  e {b|b ! T /l/ 2K  (a’ (F  -G  )+a,+(P  +G  )) 
z  XL  V  V  V  V  v  V  V 


+  btb2J/l/2Ku(ayF*+G*)+a;+(F*-G*))}  , 


where 


and 


ttr 


and 

(5.9) 


(5.11) 


(5.12) 


' 
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In  Eq.  (5.1l)  the  summation  is  over  the  "transformed" 
eigenstates  of  the  photon.  The  new  infinite  zero-point  energy 
may  be  subtracted  from  .  We  notice  that  the  first  three 

terms  in  are  very  similar  to  H0  of  (4.6)  but  of 

(4.6)  has  a  different  dependence  on  the  creation  and  annihi¬ 
lation  operators.  The  form  of  H  implies  that  the  electron 
in  level  1  or  2  can  absorb  or  emit  an  arbitrary  number 
of  "transformed"  photons  of  eigenstates  V  .  With  the  Hamil¬ 
tonian  (4.6),  the  bound  electron  could  only  emit  or  absorb 
one  "free"  photon  of  eigenstate  A.  .  This  means  that  a  "free" 
photon  corresponds  to  a.  superpoi/fron  of  "transformed"  photons,  which 
was  actually  implied  by  (5*9)*  He  should  mention  that  if 
we  treat  (5.1)  in  another  way,  i.e.  replacing  A^(r,  t )  0) , 

which  would  have  no  direct  physical  meaning  like  that  of  (5.2). 
The  H2  term  will  modify  the  energy  levels  of  the  bound 
electron  and  the  electron  operators,  b^  ^  and  t-f  2.  0  We 
can  say,  similarly,  that  a  "free"  electron,  bound  in  a 
nuclear  field  but  free  from  ,  corresponds  to  a  cloud  of 

"transformed"  electrons. 

3 .  Physical  and  Bare  Vacua  for  the  Photon 

We  shall  call  the  "transformed"  photons  physical 

* 

and  the  "free"  photons  bare.  The  new  zero-point  energy 
£f^/2,  will  be  subtracted  from  the  Hamiltonian  HT  ,  instead 
of  subtracting  £  k^/2  as  was  done  in  (4.6).  Let  us  define 


::09[r  -I"  ,  i  im  ,  H  "io  r:.l  erfl  c  ^]C 

obnoQ  as^-xoo.  t  K  njoni  93*1!  $1  <J  C 
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the  physical  vacuum  state  as  that  state  with  no 

physical  photon  but  one  electron  at  level  1  or  2.  The 
presence  of  one  electron  in  |  is  to  satisfy  the  condition 

given  by  (5«3)j  so  that  H  ^  in  (5«7)  will  not  vanish  identi¬ 
cally,  otherwise  we  would  have  the  same  Hamiltonian  as  (4.6), 
For  definiteness  we  shall  put  this  electron  at  level  2. 

Then  the  physical  vacuum  state  for  the  photon  is  defined  by 

a^|*J>  =  0  ,  (5.13) 

for  all  V  .  Similarly,  we  can  also  define  a  bare  vacuum 
state  |  §o')  for  the  photon  as 

a^  |  <i>o >  =  0  ,  (5.14) 

which  is  a  state  with  no  photon  but  one  electron  at  level  2. 

Let  us  assume  that  initially  there  is  no  perturbation 
in  our  system  so  we  can  make  up  our  initial  state  by 
applying  a^_  to  |  $0">  .  Recall,  from  Section  4  of  last 

chapter,  that  the  mean  energy  and  the  rate  of  energy  dissipa¬ 
tion  should  be  finite  for  a  physically  acceptable  state.  It 

is,  therefore,  sufficient  to  test  that  if  < §0  I  H ^  I  $0^  and 

£ 

h;  i$.>  are  finite,  our  initial  state  will  be  physically 
acceptable.  Rewriting  H*  from  (5.1l)  and  (5.8),  and  then 


on  rftflw  9ctB J  as  <o$>j  eJs  '  so*‘-*3  '-fii,i. 

, 

STS 

-  ; 
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expanding  A  by  (4.2)  in  terms  of  operators  a^_  and  a^  s 
a  straightforward  calculation  yields 


<$o  |  |  >  =  E2  +  ^kA/2-lKv/2+(e2TT/m)  Jl/kx  .  (5.15) 


The  three  sums  are  separately  divergent  at  high  frequencies, 

the  first  two  quartically,  the  last  quadratically .  We  cannot 

say  offhand  that  the  divergences  completely  compensate  each 

16 

other  at  high  frequencies  ,  since  Kv  cannot  be  solved 
simply.  However,,  there  is  still  hope  that  (5*15)  may  be 
finite  due  to  the  difference  in  signs  for  the  three  sums. 

For  the  case  of  a  point  charge  distribution,  |rp( r)|  -  <d(r), 

we  have 


H2  =  ( e2  2ir/m)  J  £,-e^,q^t  , 


A  A 


,  A  A  ’  A  A 


(5.16) 


which  can  be  diagonalized  with  HreLc|  by  Van  Kampen's  method 
Instead  of  (A.ll),  we  have  a  simple  equation  for  the  eigen- 


17 


values  , 


( e2  lir/m)  £  l/(K2-k2)  =  1  . 


(5.17) 


, 
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Since  we  want  to  examine  the  high-frequency  behaviour  of 


frequencies  (5.17)  will  be  dominated  by  the  term  + 

where  8  is  a  small  positive  frequency.  With  the  above 
approximations,  we  obtain 


K  k,  +  e  227T/(mkA) 
v  A  A 


(5.18) 


at  high  frequencies.  We  can  see  from  (5.18)  and  (5*15)  that 
the  divergences  in  (5.15)  exactly  compensate  each  other  at 
high  frequencies.  Hence,  the  mean  energy  for  the  bare  vacuum 
is  finite  if  2  K-p/2  is  subtracted  from  the  Hamiltonian, 
instead  of  subtracting  2  k^/2  as  is  done  usually.  Let  us 
proceed  to  evaluate 


<fo  [Hy  |$o>  =  (<$o  |Hpto>)2  +  e2£|Sx|2/(2kA)  + 


/  ~ ^  \  I  9  I  9 

ip(r )  I  dr  I 


(5.19) 


■ 


. 

' 


87 


The  last  term  diverges  for  point  charge  distribution,  while 

the  second  term  converges  due  to  |  (B^|^  .  We  can  say  that 

states  built  up  from  the  bare  vacuum  is  not  physically 

acceptable  at  least  for  the  point  charge  distribution.  If 

we  assume  a  spread-out  charge  distribution  concentrated  in 

2. 

a  region  of  radius  l/k  about  the  origin,  then  |J . \  in 

(5*19)  will  provide  a  cutoff  for  frequencies  >  kc  to  make 
the  last  term  finite.  However,  there  is  no  guarantee  that 
(5.15)  will  remain  finite  for  a  spread-out  charge  distribution. 
The  arguments  above  lead  us  to  conclude  that  the  concept  of 
the  bare  vacuum  is  not  physically  acceptable.  This  is 
connected  to  the  well  known  situation  in  quantum  field  theory 
that  one  cannot  separate  out  a  bare  particle  from  the  physical 
particles,  since  a  physical  particle  always  associates  with 
a  cloud  of  bare  particles.  On  the  other  hand,  we  can  show 
that,  at  t  —  0, 


«H  |H^|$i>  =  E2 


and 


<4>J  1  H^,2  |  $}>  =  (E2  )  2  +  e2£|Fv-Gv|2/(2Kv) 


(5.20) 


. 
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which  are  both  finite.  Hence,  we  should  prepare  our  initial 
states  from  the  physical  vacuum. 
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Chapter  6.  SUMMARY  AND  CONCLUSION 

We  have  studied  two  examples  of  quantum  mechanical 
dissipative  systems.  They  are  subsystems  of  conservative 
systems.  In  the  first  example,  the  conservative  system 
consists  of  two  interacting  particles  which  can  make 
transitions  to  n  different  final  states  (real  or  virtual). 
A  single  channel  in  this  case  will  be  the  (dissipative) 
subsystem.  Here,  assuming  a  nonlocal  separable  potential, 
we  can  obtain  exact  result  for  the  partial-wave  amplitude 
which  contains  the  usual  kinematical  cuts  from  the  unitarity 
condition  and  dynamical  cuts  from  the  interaction.  However, 
there  is  no  simple  way  to  generalize  the  one-channel  N/D 
method  to  the  multichannel  problem.  Of  course,  in  our 
solvable  example,  there  is  a  matrix  A  which  has  all  the 
properties  of  the  desirable  D  matrix,  but  the  complete 

expression  for  the  partial-wave  amplitude  cannot  be  written 

-1 

as  a  quotient  of  two  matrices,  ND  .  Bound  states  and 
resonances  are  given  by  zeros  of  det(A)  as  expected. 

The  second  example  concerns  a  solvable  model  in 
nonrelativist ic  electrodynamics.  Results  of  level  shift 
and  line  breadth,  corresponding  to  a  damped  classical 
point  charge,  are  obtained.  The  decay  behaviour  of  a  state 
has  the  usual  exponential  form  for  small  time  (^  half  life 


1  «.  -1  ; 


■ 
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of  the  state)  .  This  damped  exponential  and  the  line  breadth 
can  be  considered  typical  characteristics  of  a  quantum 
mechanical  dissipative  system. 

In  the  last  chapter,  we  have  considered  modification 

of  the  second  example,  which  changes  a  bare  photon  to  a  superposition 

of  physical  photons  as  implied  by  (5.9).  We  have  also  argued 

that,  for  physically  acceptable  states,  the  transformed  zero- 

point  energy  Jk^/2  should  be  subtracted  rather  than  the 

radiation  field  (free)  zero-point  energy  Jk^/2.  One  notices 

that  IKV/2  arises  from  the  Heisenberg's  equations  of 

motion  for  a'  ,  i.e. 

v  5 


ayt)  =  i[HT,  a^(t)]  ,  (6.1) 

by  the  usual  commutation  relation 

[a^(t) ,  a^l(t)]  =  5vv,  .  (6.2) 

On  the  other  hand,  the  equations  of  motion  for  the  bare  pho¬ 
ton  operators , 


ax(t)  =  i[HT,  ax(t)] 


(6.3) 


is  obtained  by  assuming 


. 


■ 
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[ ( t ) 3  aj , ( t ) ]  =  6aa,  .  (6.4) 

An  entirely  different  zero-point  energy  £k^/2  arises. 

If  we  demand  that  (6.3)  should  also  give  the  physically 
acceptable  result,  namely  the  zero-point  energy  £k^/2  , 
then  the  commutation  relation  (6.4)  will  no  longer  hold 
and  it  should  be  modified.  This  leads  us  back  to  Wigner’s 
result  that  from  Heisenberg's  equations  of  motion  (6.3), 
the  zero-point  energy  Jk^/2  can  be  obtained  if  one  assumes 
a  different  commutation  relation. 
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APPENDIX  DIAGONALIZATION  OF  Hrad  +  H2 

In  this  appendix  we  shall  deal  with  the  diagonaliza- 
tion  of  Hrad +  ^2  in  some  detail.  Substitution  of  (5.4)  into 
(5-7)  yields 

Hrad  +  H*  =  (1/27(PA+kAqf  +  (eV2  )  7  ,  (gA  A  ’ p  XPX  '  +  fXX  'qXqX  •  + 

A  A  A 

+  hA,APxqA,+hXA,qAPA,)  , 

(  A  .  1 ) 


where 


r 


■XX 


(  Air/mk^k^  ,  *e^  t 


sin(  k,  *r )  sin(ic,  ,  •  r )  |  ip  |  2dr  = 


X 


‘X  ’X 


r 


fXXT 


(  4tt/iti) e^  •  e^  , 


cos  ( k^  •r )  cos  ( k^  ,  *r )  |  ijj  |  2dr  = 


f 

X  'X  9 


hXX' 


=  (-Hir/mk^  ,  )ex^x, 


cos(k^*r)sin(k^  f  *r)  1 |  2 dr 


(A. 2) 


, 


' 
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We  shall  consider  (A.l)  as  a  classical  Hamiltonian  of  finite 
degrees  of  freedom,  i.e.A.=  l,2,  ...  N.  Assume  that  the  final 
result  can  be  carried  through  for  if  all  sums  remain 

finite . 


First  let  us  eliminate  the  cross  terms  in  p.  and  q. 

A-  A. 

by  the  following  transformation: 


P  A  ^A,(^AA,PA’  +  ^AA,qA,) 

qA  =  ^A'^AA'Pa'  +  °^AA  ’  qA  '  ^  *  (A. 3) 


The  transformation  matrices,  ,6,6,  and JD  ,  are  deter- 
mined  by  the  two  requirements: 

(a)  transformation  (A. 3)  is  canonical,  viz.  the  fundamental 

Poisson  bracket  relations  are  invariant  under  canonical 

1 

transformations;  and 

(b)  cross  terms  in  (A.l)  vanish  after  the  transformation. 
Condition  (a)  leads  to 

£a"®AA"^VA"  “  <^AA"^A,An)  =  6AA' 

A\"  “  ^AA"^A'An)  =  0 


. 

. 


. 

. 
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-  JO 


XX 


^Vxm) 


0, 


(A. 4) 


while  the  second  condition  demands 


h(^AA2^AA2  +  kAC’AA1^AA2)  +  ehAA  .  (Sxx  •  ^AA  ,  S\  '  A2  + 


+  fAA'  <^AAi^)A,A2+hA,A^AAi"^A,A2+hAA'^AAi1^A,A2^  0 


(A. 5) 


Equations  (A. 4)  and  (A. 5),  in  principle,  can  be  solved  and 
thus  determine  transformation  (A. 3) .  Then  (A.l)  becomes 


H 


rad 


+  H; 


=  ^a.aAA 


a  Pa  Pa  A  P  A  Pa  ) 

1  A2  A  1  A  2  A  1  A  2  A  1  A2 


(A. 6) 


where 


^XiX2  ^X^XXi^XX2+kX  ^  XXi^XX^  +  + 


XX,V&XX'  XXi  X’X; 


+  fXX’^XX1^X'X2  +  hX'X^XX1^X*X2  +  hXX'^  XXl^X'\2)  ~^X2X1 


. 


. 
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and 


^AiA2  \\fi  \\2  +  ^  A  ^A  A  i  ^  A  A  2  ^  +  6  ^AA  '  (gA  A  T  ^  A  A  i  ^  A  '  A  2  + 


+  ^A  A  AA  A  ’  A2  +  hA’ A^  AAi^A’A2  +  hAA  '  ^  AA  ,  ^A  '  A?  }  ^A2Ai  * 


(A. 7) 


The  sum  of  ^ra(j  +  H2  P°sit ive-def init e  and  real, 

as  should  be  expected  from  the  radiation  energy  plus  a  per¬ 
turbation  energy  (  (5*7)  and  (A.l)  ).  Then  (A. 6)  consists  of 

two  real  positive-definite  quadratic  forms  which  can  be  dia- 

48 


gonalized  simultaneously 


The  canonical  transformation 


is 


t  — 


A 


=  It 


v 


,  p 
Av  v 


‘A 


=  I(T  b 


V 


vAQv 


(A. 8) 


which  satisfies  condition  (a).  In  order  that  (A. 8) 
(A. 6)  to  a  form  like  Hrad  3  the  first  sum  in  (A.l) 


reduces 
we  require 


and 


. 


■ 


b 
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WT  Kx^XX'^  \'X'  Kv6vv’  J  (A. 9) 

where  are  the  eigenvalues  to  be  determined.  From  (A. 9), 

we  find 

h'(%A'  -  Kv(A."1)u,)(T-1)vA,  =  0  ,  (A. 10) 

which  may  be  considered  as  a  set  of  N  homogeneous  algebraic 

equations  for  the  transformation  matrix  elements  (T_l  )  , 

VA 

where  V  is  fixed.  The  necessary  and  sufficient  condition 

that  these  equations  have  a  solution  is  that  the  determinant 

z 

of  their  coefficients  vanishes,  i.e.  the  eigenvalues 
are  the  roots  of  the  secular  equation 


det 


^i:l-k2(A  b11,  th12-k2(X,  1)12,  •••  b 

]?t21-K2(<C"1)21,  )TU,2-k2(A'1)22J"  •  ^n-k2(X-_1) 


IN 


2N 


^N1  K  ^Nl5 


^N2  K  ^ N2  5  *  '  %fN  K  ^ 


NN 


=  0 


(A. 11) 


ra '!  j  ' 

, 
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As  N-^00  ,  we  have  a  secular  determinant  of  infinite  dimension 
in  (A.ll).  To  obtain  such  a  determinant,  mathematical  induc¬ 
tion  is  employed  to  determine  a  general  solution  for  arbitrary 

4  9 

dimension  N  and  then  N  is  allowed  to  approach  infinity. 

47 

Another  method  to  diagonalize  (A. 6)  is  given  by  Van  Kampen 
for  the  special  case  that  ^  and  ~  MA^AA.' +  r\rrY  ‘ 

We  obtain  finally,  from  (A. 6),  (A. 8),  and  (A.9)j 


H 


rad 


N 

ib 

V  =  1 


+  H2  =  3§J(  P2  +  K2Q2)  , 
L  v  V  V 


(A. 12) 


z  .  . 

where  K_^>  ,  as  determined  from  (A.ll),  are  positive  and  real 

since  we  have  assumed  that  Hrad ^  ^2  Posikive  and  real. 

2. 

For  simplicity,  we  shall  also  assume  that  are  all  distinct 

In  (A. 12)  we  can  let  N-»oo  as  the  diagonalized  form  of  Hra(j  +  . 

Quantization  can  be  carried  out  to  yield 


had  +  H*  “  IKvK+V  +  . 

v 


(A. 13) 


where  a^"  and  a^  are  the  new  set  of  creation  and  annihi¬ 
lation  operators  defined  similar  to  (5*6).  A  new  infinite 
zero-point  energy,  ^  K^/2,  appears.  One  should  note  that 


K2  =  k2  for  e2  =  0  , 
v  A  s 


(A. 14) 


as  one  compares  (A. 12)  with  (A.l). 
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